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ABSTRACT 


The  current  analytical  models  of  impulsive  noise  and 
system  response  thereto  are  found  numerically  intractable, 
indicating  a  need  for  engineering  approximations.  One,  based 
on  finding  an  "equivalent"  unimodal  system  having  approximately 
the  same  response  s tatis tics--here  the  crest  statistics-- 
is  examined,  it  is  found  that  fair  agreement  obtains,  making 
it  possible  to  use  two  reasonably  accessible  parameters 
(v  and  ui  )  to  characterize  any  linear  multimodal  system 
response  and  thus  predict  its  crest  statistics  systematically, 
provided  that  "standard"  curves  of  unimodal  response  crest 
statistics  are  available  for  statistically  identical 
forcing.  Limitations  of  this  promising  approach  are  explored, 
encouraging  its  use  in  reliability  prediction  pending 
further  studies.  Also,  the  need  for  further  studies  of 
impulsive  noise  response  crest  distribution  "tails"  is 
noted  for  prediction  of  overload  and  wearout  reliability 
on  the  basis  of  wear-dependent  failure  rate  concepts,  since 
extrapolation  of  empirical  data  into  the  "tail"  (here  done 
on  the  basis  of  a  Weibull  crest  fit)  is  increasingly  necessary 
the  more  "impulsive"  the  response  is  (i.e.,  the  more  the 
parameter  v  is  less  than  unity)  . 
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INTRODUCTION 


The  prediction  of  system  reliability  under  threat  of 
overload  and  wearout  failure  is  currently  performed  in  term* 
of  system  response  "peak*  (crest,  maxima,  or  rise) 
statistics  for  the  standard  models  of  these  two  basic 
failure  mechanisms.  Hence,  one  important  problem  confronting 
the  systems  analyst  is  to  find  these  "peak"  statistics — 
particularly  the  level  crossing  (crest)  statistics  which 
yield  conservative  estimates--for  the  system  in  question 
under  an  anticipated  random  load. 

Due  to  the  prevalence  of  nongaussian,  particularly 
impulsive  (characterized  by  large-amplitude  bursts  against 
a  low-amplitude  background)  noise  in  •'^alistic  models  of 
random  load  environments,  the  analytic  ;.x  difficulties  are 
formidable  in  an  accurate  computation  of  the  pertinent 
"peak"  statistical  distributions.  Hence,  one  would  like 
to  learn  how  to  apply  current  analytical  models  for 
statistical  system  response  to  impulsive  noise  in  order  to 
estimate  "peak"  statistics.  Moreover,  one  desires 
"engineering"  response  models  which  will  serve,  in  particular 
physical  situations,  until  more  general  analytical  results 
have  been  found.  Thus,  it  is  of  more  than  academic  interest 
to  learn  both  the  limitations  and  potentialities  of  the 
results  obtained  for  the  impulsive-noise-response  of  a 
simple  linear  harmonic  oscillator,  mutatis  mutandis,  to 
predict  the  “peak"  statistics  (and  thus  reliability)  of  more 
complex  multimodal  linear  systems. 

This  report  will  first  briefly  summarize  the  state  of 
the  art  in  the  analysis  of  impulsive  noise  response  of  linear, 
time- invariant  systems.  Then,  possibilities  of  finding 
approximately  "equivalent"  unimodal  systems — facilitating 
reliability  prediction  for  complex  multimodal  systems — 
will  be  explored.  It  will  be  seen  that  multimodal  system 
"peak"  statistics  may  often  be  conservatively  (but  not 
excessively  so)  estimated  for  reliability  predictions,  the 
system  being  characterized  by  two  physically  meaningful 
paramecers,  a  "time  constant"  t  (or  bandwidth)  and  an 
"equivalent  resonant  frequency”  uu  . 
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It  was  found  that  utilizing  initial-slope  decay  times 
vielded  lower  bound  estimates  for  the  empirical  crest 
distributions,  while  the  bandwidth  parameter  yielded 
conservative  (upper  bound)  estimates.  Predictions  of  the 
exact  frequency  of  zero-crossinns  (necessary  for  normal¬ 
ization)  ,  as  well  as  the  eauivalent  unimodal  system  resonant 
frenuenev,  presented  certain  difficulties;  but,  one  obtains 
normalized  crest  distribution  predictions  (frequency  of 
level-crossings  divided  by  that  of  zero-crossings)  that  are 
above  the  corresponding  unimodal  statistics  bv  a  factor  of 
the  order  of  unitv.  Potentialities  and  limitations  of  these 
models  are  elucidated,  and  areas  of  further  research  delineated 
in  suhseauent  sections  of  this  report. 
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II 


STATISTICAL  MODELS  FOR  SOURCES  AND  SYSTEMS 


Reliability  predictions  for  a  complex  system  under  random 
loading  generally  assume  knowledge  of  a  relation  between 
some  particular  aspect  of  the  system's  response  and  system 
failure  (ref.  1) .  With  this  information  one  then  usually 
attempts  to  predict  the  average  time  to  failure.  For  this 
purpose  the  system  response  "peak"  distributions  (crests, 
maxima,  or  rises)  are  typically  utilized  in  overload  and 
wearout  reliability  calculations.  Previous  workers  have 
found  that,  insofar  as  partial  wear  (damage)  predictions 
are  concerned,  a  hierarchy  exists  among  the  "peak" 
statistics  (refs.  2,  3)  .  For  a  given  situation  the  crest 
statistics  offer  the  most  conservative  (i.e.,  shortest) 
estimate  of  the  time  to  failure,  while  the  rise  statistics 
provide  the  most  optimistic  predictions.  For  the  purposes 
of  this  report  the  "conservative"  crest  statistics  will  be 
the  main  statistical  distribution  of  interest.  The  crest 
statistics,  of  course,  are  also  of  particular  interest  in 
overload  failure  prediction  (refs.  4,  5). 

The  response  "peak"  statistics  may  be  approached  from 
either  analytical  or  empirical  viewpoints.  For  analytical 
predictions  one  must  have  sufficient  "higher  order" 
statistical  information  (cf.  Section  III)  about  the  noise 
source  as  well  as  a  characterization  of  the  system  under 
study  so  that  formal  expressions  for  the  system  response 
statistics  can  be  obtained.  This  approach,  if  successful, 
permits  a  systematic  study  of  alterations  in  the  average  time 
to  system  failure  precipitated  by  changes  in  either  the 
noise  source  parameters  or  the  system  parameters.  However, 
in  most  physical  situations,  particularly  nongaussian  noise 
loading,  when  such  detailed  information  is  not  available,  one 
must  resort  to  fitting  empirically  the  system  response  "peak" 
distributions  by  analytical  curves  (say  Weibull  (ref.  4,  5), 
exponential,  Gamma,  etc.).  At  best,  one  might  observe  the 
development  of  trends  which  could  be  linked  directly  to 
parameter  changes  either  in  the  noise  source  or  in  the  system 
itself.  For  the  sake  of  general  applicability  the  analytical 
approach  outlined  above  will  be  attempted  as  extensively  as 
possible.  However,  when  an  impasse  is  reached,  the  second 
method  will  be  utilized  without  further  justification. 
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1. 


IMPULSIVE  NOISE  MODELS 


In  an  attempt  to  analytically  characterize  environmental 
impulse  noise,  various  models  have  been  proposed  and  have 
received  study  in  the  literature  (refs.  6-9) .  Generally 
speaking,  these  models  possess  a  two- fold  purpose.  First, 
they  strive  to  describe  as  accurately  as  possible  the 
environmental  noise,  implicit  in  this  operation  is  the 
establishment  or  the  definition  of  the  parameters  which  are 
deemed  crucial  in  describing  the  physical  noise  process. 
Second,  the  models  are  formulated  so  as  to  facilitate 
theoretical  predictions  of  the  desired  noise  statistics. 
Presumably,  then,  the  final  analytical  results  (e.g.,  system 
response  crest  or  maxima  statistics  or  the  response  first- 
order  probability  density)  are  explicitly  dependent  on  the 
parameters  chosen  to  describe  the  noise  process. 

When  analyzing  a  system  forced  by  impulse  noise  one  has 
two  alternatives.  Either  the  actual  noise  source  or  the 
response  of  the  system  itself  stay  be  related  to  an  impulse 
noise  model.  With  the  first  approach,  i.e.,  noise  source 
modeling,  one  is  still  concerned  with  calculating  the 
response  of  sosse  system  coupled  to  that  source.  Thus,  it 
is  also  necessary  that  the  system  be  characterized. 

Typically  what  is  needed  is  some  “input- output"  transfer 
information  of  the  system  to  account  for  its  response  to  the 
given  noise  source.  Although  this  approach  may  lead  to  more 
complex  analytical  operations,  it  does  subdivide  the 
pertinent  environmental  parameters  into  two  groups,  one 
associated  with  the  noise  source  and  the  other  group  related 
to  the  system  itself.  In  this  fashion  the  system  response 
can  be  directly  related  to  changes  in  the  noise  source  or 
to  changes  in  the  system.  The  other  alternative,  that  of 
modeling  the  actual  system  response  directly,  may  lead  to 
easier  analytical  analysis.  But  at  the  same  time,  one  loses 
the  ability  to  isolate  any  "cause-effect"  relationships  which 
may  exist  between  the  response  statistics  and  the  individual 
parameters  of  the  noise  source  and  of  the  system. 

Two  analytical  models  are  used  below.  The  first,  the 
Poisson  impulse  noise  model,  is  used  here  to  model  the  actual 
noise  source  (and,  mutatis  mutandis,  the  response) .  The 
second  approach.  Hall's  generalized- t  model,  must  be  applied 
directly  to  the  system  response  under  impulse  loading. 
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a.  Poisson  Impulsive  Noise  Model 


S.  0.  Rice's  Poisson  shot  noise  model  (refs.  6,  8,  9) 
offers  the  greatest  promise  insofar  as  analytical  elegance 
is  concerned  while  at  the  same  time  elucidating  the  physically 
meaningful  parameters.  Mathematically,  this  model  represents 
the  forcing  noise  y(t)  as  a  "train"  of  impulse  functions: 

y ( t)  *  ai  6(t-ti)  .  (1) 

i 

The  occurrence  times  -1 1 .  j  are  assumed  i)  to  be  statistically 
independent  of  the  impulse  strengths  la.),  and  ii)  to  occur 
in  a  Poisson  fashion*  with  averago  impulse  rate  (or  frequency) 
v.  The  a .  ,  in  turn,  are  assumed  to  be  mutually  independent 
random  variables  with  a  common  impulse  strength  probability 
density  pa(a).  Unfortunately,  the  resultant  formal  quadrature 
or  series  expression  for  say  the  response  first-order 
probability  density  or  crest  statistics  of  a  linear  system 
forced  by  Poisson  noise  are  usually  not  readily  converted  into 
numerical  results  (refs.  ID,  11)  suitable  for  reliability 
predictions.  Typically,  formal  results  contain  either 
characteristic  functions  which  cannot  be  readily  inverse- 
transformed  or  infinite  summations  of  Hermite  polynomials 
whose  convergence  is  either  extremely  slow  or  even  entirely 
in  doubt.  Thus,  utilization  of  the  Poisson  model  which  is 
physically  appealing,  ia  somewhat  limited  pending  advances 
in  numerical  computational  techniques. 

b.  Generalized- t  Impulsive  Noise  Model 

The  genera lized-t  model  reported  in  a  recent  work 
(ref.  7)  by  H.  M.  Hall  also  offers  promise  for  characterization 
of  impulse  noise.  Actually,  Hall  was  primarily  concerned 
with  modeling  radio  noise  caused  by  atmospheric  disturbances 
(e.g.,  lightning)  as  observed  at  the  output  of  a  communications 
receiver.  Thus,  his  results  are  most  directly  applicable  to 
filtered  "impulsive"  atmospheric  radio  noise.  Hall's  model 


The  Poisson  distribution  giving  the  probability  of  n  impulses 
occurring  in  an  observation  time  T  is  given  by 


PT<n) 


( vT) n  @~vT 


is  therefore  of  the  second  type  outlined  above,  i.e.,  one 
associated  directly  with  the  system  response  process.  From 
a  mathematical  point  of  view  Hall's  model  is  applicable 
without  major  modification  only  to  the  output  response  of  the 
system  being  forced  by  impulsive  noise. 


A  formulation  of  the  generalized-t  model  proceeds  as 
follows.  Let  the  response  process  y(t)  be  represented  as 
a  product  of  two  statistically  independent  noise  processes 
denoted  by  m(t)  and  n(t)  ,  i.e., 

y ( t)  »  m(t)n(t)  .  (2) 


Hall's  analysis  required  that  the  "carrier"  n(t)  be  an 
ensemble  member  of  a  narrowband  Guassian  process  with 
correlation  function  Rn(t) »  while  the  "modulating"  process 
m(t)  be  a  slowly  varying  (with  respect  to  n(t))  process 
with  a  known  first-order  probability  density  function,  p^m) . 
Actually,  n(t)  need  not  be  restricted  to  a  narrow-band 
process.  It  turns  out  that  a  necessary  requirement  is  that 
the  frequency  spectra  of  n(t)  and  m(t)  be  non-overlapping. 


Using  the  above  assumptions  about  m(t)  and  n(t)  one 
can  formally  express  various  densities  of  y(t)  as  integrals 
(i.e.,  averages)  of  known  Gaussian  densities.  For  example, 
one  immediate  result2  is  that  the  first-order  probability 
density  of  y(t),  py(y) ,  be  given  by 

oo 

P  (y)  ■  f  T^T  P  Ap  (x)dx  •  (3) 

ry  1  J  |x|  rn  x  *m 

-oo 


Since  n(t)  is  a  member  of  a  Gaussian  process,  the  form  of 
Pjjt*)  is  known.  Thus,  specif i  :at ion  of  Pn(")  formally  defines 
Py(y) •  Of  course,  analytical  evaluation  of  the  above  integral 
may  become  an  obstacle.  Other  statistical  distributions,  for 
example  the  crest  or  maximum  statistics  of  y(t),  also  can 
be  obtained,  at  least  in  quadrature;  such  problems  will  receive 
further  attention  in  Section  HI.  Thus,  the  formal  analytical 
simplicity  of  Hall's  impulse  noise  model  is  appealing. 


2 

The  requirement  that  m(t)  be  "slowly  varying’’  is  not  necessary 
in  order  to  obtain  this  particular  result  for  Py(y),  but  this 
restriction  is  quite  important  in  computing  "peak"  statistics. 
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although  its  suitability  for  engineering  predictions  is  still 
a  point  in  question. 

Unfortunately,  a  logical  association  between  the  pertinent 
parameters  of  Hall's  model  and  available  environmental  data 
is  quite  obscure.  In  contrast,  one  of  the  attributes  of  the 
Poisson  model  is  the  rather  obvious  correspondence  between 
the  noise  model  impulse  rate  v  and  the  observed  rate  of  occur¬ 
rence  of  noise  pulses  from  field  data  (ref.  12) .  In  his 
work.  Hall  modeled  response  envelope  data  obtained  by 
filtering  atmospheric  radio  noise.  He  found  that  a  chi- 
squared  distribution  for  the  random  variable  b  fc  1/m, 

u  M/2 
(— ) 

2'  i  M- 1  M  .  2i 

Pb<b)  *  tt  lbi  expr  — J  b  }  *  (4) 

b  o"r(|)  1  2  c?  J 

yields  a  proper  "asymptotic  behavior"  of  the  first-order 
probability  density  py(y)  for  large  y  as  suggested  by  other 
investigators  (ref.  13) .  Then,  by  adjusting  the  two  para¬ 
meters  in  Eq.  (4)  namely  (M  and  o) ,  one  is  able  to  obtain 
"best  fits"  to  filtered  atmospheric  radio  noise  envelope 
data  (first-order  envelope  probability  densities  and  envelope 
level  crossing  rates) .  The  selection  of  appropriate  model 
parameters  is,  then,  based  solely  on  data  observed  at  the 
output  of  the  sys tern  in  question.  Thus,  although  one  may 
be  able  to  model  the  system  response,  it  is  not  at  all  obvious 
how  these  analytical  response  statistics  are  related  to 
certain  parameters  of  either  the  noise  source  (e.g.,  average 
impulse  occurrence  rate)  or  of  the  system  itself  (e.g., 
the  system  bandwidth) . 

One  faces  comparable  problems  when  endeavoring  to  obtain 
damage  estimates  for  any  physical  situation  when  employing 
the  generalized-t  model  to  characterize  system  response  to 
impulsive  noise  environments.  Moreover,  the  use  of  this 
model  in  damage  calculations  is  misleading  in  that  the 
formal  results  may  produce  infinite  damage  moments  (cf. 

Eq.  (30)).  Nevertheless,  on  consideration  of  physically 
useful  models  for  impulsive  noise,  the  above  mentioned 
limitations  of  the  generalized-t  model  certainly  are  not 
so  formidable  as  to  rule  out  the  method,  mutatis  mutandis. 


2. 


SYSTEM  MODELING 


Whan  using  the  Poisson  model  of  the  preceding  section, 
one  must  necessarily  characterize  analytically  the  linear 
sys  tem  under  study.  Even  though  the  generalized-t  model 
is  applicable  only  to  the  system  response,  one  would  suspect 
that  a  detailed  characterization  of  the  linear  system  might 
possibly  be  directly  associated  with  resultant  changes  in 
the  response  statistics.  Thus,  even  with  Hall's  modeling 
procedure,  system  modeling  should  prove  quite  useful. 

None  of  the  preceding  discussions  have  focused  directly 
on  the  form  of  system  characterization  to  be  used.  The 
formulation  of  the  Poisson  noise  model  in  terms  of  "impulse 
functions"  suggests  strongly  that  the  system  impulse  response 
function  h(T)  should  be  embodied  in  the  system  modeling 
process.  Once  a  h (t)  is  specified,  the  system  in  question 
(assumed  in  this  report  to  be  linear  and  time-invariant) 
is  fully  described  (ref.  14) .  For  a  complex  system,  however, 
one  may  not  be  able  to  obtain  h(t)  in  a  tractable  form.  In 
fact  only  for  rather  simple  systems  can  h(i)  be  readily 
expressed  in  closed  form. 


For  our  purposes  perhaps  the  most  common  system  for 
which  this  characterization  is  possible  is  the  simple  linear 
harmonic  oscillator.  It  is  a  natural  building  block  for 
more  coi^plex  linear  systems:  of  particular  interest  here  are 
underdamped  (high  Q)  systems.  Such  a  system  has  an  impulse 
response  (ref.  15)  given  by 

<x> 

o 


h  (T) 


A  e 


20 


T 

sinu)  x 
o 


(5) 


where  oo0  is  the  undamped  resonant  frequency,  Q  is  the  system 
quality  factor,  and  A  is  a  normalization  ("gain")  factor, 
in  this  form,  Eq.  (5)  is  particularly  amenable  to  analytical 
calculations  and  predictions  of  system  reliability. 

Characterization  of  more  complex  systems  by  the  speci¬ 
fication  of  h(T)  is  not  readily  accomplished.  For  example, 
a  two-degree-of- freedom  system  which  can  be  described  by 
a  pair  of  coupled  second-order  linear  differential  equations 
(with  constant  coefficients)  such  as  Eq.  (46)  in  Section  V 
yields  an  impulse  response  of  the  form 


8 


-v 

h(T)  «  A^t  Sinju^l)  + 


"V 

8in(^2T) 


(6) 


for  certain  choices  of  system  parameters.  It  should  be 
observed  that  Eq.  (6)  basically  does  involve  linear  combina¬ 
tions  of  the  more  simple  h(t)  characteristics  of  the  single¬ 
mode  case,  Eq.  (b)  .  However,  the  analytical  form  of  Eq.  (6) 
is  still  sufficiently  complex  to  obviate  convenient 
utilization  for  predicting  system  response  statistics  to 
impulsive  loading.  If  one  considers  higher-degree-of- 
freedom  systems,  the  situation  is  further  obscured  and 
for  a  distributed- parameter  system,  the  expressions  for 
h(T)  are  available  only  as  infinite  sums  over  the  tiigenmodes 
of  the  system  (cf.  discussion  of  a  clamped  circular  plate 
in  the  Appendix) . 

Faced  with  these  inherent  difficulties  one  is  forced 
to  investigate  as  to  whether  or  not  less  information  than 
the  complete  specification  of  h(f)  would  suffice  for  making 
damage  predictions.  As  an  example,  one  might  ask  if  some 
"average  system  response  time"  information  is  sufficient; 
or  possibly  an  "equivalent  aingle-degree-of-freedom  system" 
could  be  found  and  used  as  a  statistically  reliable  model  for 
the  more  complex  system  when  making  the  predictions.  These 
specific  questions  and  others  will  be  discussed  later  in  more 
detail.  It  will  be  shown  that  for  certain  engineering 
estimates  one  can  extend  known  unimodal  system  behavior  to 
new  situations  involving  complex  systems  utilizing  only  a 
minimal  aincunt  of  information  about  the  system  in  question. 
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III 


"EXACT  THEORIES"  FOR  PREDICTING  "PEAK"  DISTRIBUTIONS 


The  response  "peak"  distributions  under  consideration 
formally  can  be  expressed  in  terms  of  higher-order  response 
joint  distributions.  More  often  than  not,  however,  calcula¬ 
tions  based  on  these  formal  solutions  cannot  be  converted 
into  numerical  results,  in  some  cases  computational  dif¬ 
ficulties  can  bt  circumvented  by  using  approximations  valid 
in  the  low- impulse  frequency  limit  as  discussed  in  sub-section 
1.  In  sub-section  2  "exact"  results  are  derived  for  the 
crest  and  maxima  statistics  employing  the  generalized-t 
response  model.  Although  the  analytical  results  are  relatively 
simple  in  form, their  utilization  in  damage  prediction  i3 
still  a  point  in  question;  no  relationships  between  the 
model  parameters  and  field  parameters  have  yet  been  established, 

1.  POISSON  IMPULSIVE  NOISE— LOW  FREQUENCY  LIMITS 

As  stated  earlier  formal  expressions  for  the  system 
response  crest,  maxima  and  rise  distributions  are  known 
(refs.  6,  16).  Results  for  the  crest  and  maxima  statistics 
follow  from  a  route  outlined  by  S.  O.  Rice  (ref.  6),  namely, 
that  for  a  stationary  random  process  y(t)  the  positive- slope 
level  crossing  rate  of  a  fixed  level  y  .  v_+(y  ),  is  given 


by 


oo 

vc+(yo)  “  JP  y  Py,y(yo'y)dy  ' 


(7) 


while  the  maxima  frequency  pre lability  density  m  (y  )  can  be 
expressed  as 


m+(yQ) 


o 

r 

j 

-oo 


py,y,y(yo'°'y)dy 


(8) 


( •) 


and  py,y.y( 


In  the  expressions  Eqs .  (7)  and  (8),  py 
represent  the  simultaneous  joint  probability  densities  of  the 
response  with  its  first  derivative  and  of  the  response  with 
its  first  and  second  derivative,  respectively.  The  formal 
expression  for  the  "rise"  distribution  (refs.  1,  16)  involves 
the  "method  of  inclusion  and  exclusion"  leading  to  an 
untractable  series  representation. 


.) 
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Formally,  the  crest  and  maximum  distributions  can  be 
expressed  as  integrals  of  doubly- inf inite  and  tr iply-inf inite 
series  of  generalized  Hermite  polynomials  (ref.  17)  with 
the  expansion  coefficients  related  in  known  form  (for 
Poiason  impulsive  noise  response)  to  the  noise  source 
F  rameters  (v,pa(a))  and  also  to  the  system  impulse  response 
h(x)  (ref.  18).  Unfortunately,  earlier  work  has  indicated 
(refs.  10,  11)  that  for  the  cases  of  interest,  even  a  singly- 
inf inite  Hermite  expansion  (for  the  response  first-order 
probability  density)  does  not  converge  rapidly  enough  to 
allow  numerical  computational  techniques  to  produce  usable 
results.  One  would  expect  even  greater  difficulty  with  the 
more  complex  expansions  arising  from  Eq.  (7)  and  Eq .  (8), 
even  after  the  integration  has  been  formally  performed. 
Meaningful  results  for  the  "peak"  statistics  can  be  obtained 
only  for  the  highly  nongaussian,  "low  v"  limiting  cases  which 
are  of  considerable  interest.  Two  such  approaches  are 
illustrated  below. 

The  first  approach  has  been  reported  earlier  (refs.  10, 

11,  19)  and  will  be  described  here  very  briefly.  Tf  one 

assumes  Poisson  impulsive  noise  excitation  of  a  linear  single- 
degree-of- freedom  (1DOF)  system  with  a  fairly  high  quality 
factor  (say  Q  £  5.0),  and  if  the  noise  source  impulse  rate 
v  is  sufficiently  small  (v  &  ^q/Q) .  then  the  normalized 
crest  statistics  are  given  by  (refs.  10,  19) 

oo 

c+(t|)  »  4y  [  da  p  (a)  log  ^*7^  .  (9) 

a  V  y 

2 

Here,  r|  is  the  normalized  (with  respect  to  the  rms  output) 
response  variable,  pa(a)  is  the  noise  impulse  strength 
density,  and  y  is  the  in.pulse-noise-density  parameter  defined 
by 


o 


Comparison  with  experiment  has  shown  (ref.  11)  that  Eq .  (9) 
is  valid  for  y  <  0.5.  This  prediction,  Eq.  (9),  for  the 
crest  statistics  is  valid  for  high-Q  unimodal  systems.  One 
would  like  to  extend  this  approach,  based  on  observations  of 
the  system's  impulse  response  (here  a  single  damped-s inusoid) 
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in  order  to  obtain  appropriate  limiting  forms  of  the  response 
crest  and  maximum  distributions  for  more  complex  systems. 

In  principle,  this  approach  to  predicting  response 
peak  statistics  in  the  low  impulse  frequency  limit  is 
applicable  to  all  time- invariant  linear  systems  and  (inter  alia) 
to  all  three  types  of  "peak"  statistics;  the  crests,  maxima 
and  the  rises  (ref.  3) .  It  calls  for  detailed  knowledge  of 
certain  aspects  of  the  system's  impulse  response  function 
h(r).  The  response  of  a  time- invar iant  linear  system  to 
the  Poisson  noise  of  Eq .  (1)  is: 


y(t)  «  ai’n(t-ti) 
i 


(11) 


Assuming  that  h(t)  decays  essentially  to  zero  before  some 
time  T,  and  requiring  v  «  1/T  Eq .  (11)  becomes  a  sequence 
of  essentially  nonoverlapping  impulse  responses  with 
independent,  random  scale  factors  a^.  Thus,  the  average 
exceedance  frequency  of  peaks  of  type  p  (p  =  crests,  maxima 
or  rises)  above  the  level  yQ,  vp(y0) -  is  merely  the  frequency 
of  these  peaks  in  a  single  impulse  response  with  the  scale 
factor,  a,  appropriately  averaged  over  its  distribution,  i.e., 


v  (y  ) 
p  o 


n° 


p  (a) da  , 

3 


(12) 


where  N  (yQ,a)  denotes  the  integral  number  of  peaks  above  the 
level  yQ  in  an  impulse  response  of  scale  factor  a.  The  peak 
exceedance  probability  ("peaks"  being  by  convention  positive) 
is  just  Eq.  (12)  normalized  to  unity  at  yQ  =  0,  Vp(yp) /vp(0) . 
The  only  task  remaining  then  is  the  detailed  evaluation  of 
Np(y0.  a)  . 


The  crest  distribution,  being  of  prime  concern  here,  will 
be  considered  in  detail.  Figure  1  illustrates  the  impulse 
response  of  strength  \  >  0  for  a  typical  complex  oscillatory 
system.  It  can  be  seen  from  this  figure  that  for  y0  in  any 
of  the  ranges  aM^+i  <  yQ  <  aM^  (or,  equivalently,  for 

<  a  <  y</Mi+i)  the  number  of  pos it ive- slope  crossings 
of  the  level  yQ,  Nc(y0,a),  is  constant  and  equal  to  the 
difference  between  the  number  of  maxima  and  the  number  of 
minima  above  yc;  similar  arguments  apply  for  a  ;  0.  Thus, 

Eq.  (12)  can  be  decomposed  into  a  sum  of  integrals  over  the 
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ranges  of  "a"  for  which  Nc(y0,a)  is  constant, 
the  crest  distribution  by  c+(yQ)  we  have: 


c+(y0) 


vc(0) 


N  . 
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/  M,. 


p  (a)  da 
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(13) 


In  this  expression  is  the  number  of  maxima  minuB  the 
number  of  minima  in  h(r)  greater  than  or  equal  to  for 

i  >  0  and  less  than  or  equal  to  Mi  for  i  <  0.  Also,  the 
upper  sign  in  the  t  corresponds  to  i  >  0  and  the  lower  sign 
to  i  <0.  Breaking  each  integral  up  into  the  sum  of  two 
semi- inf inite  integrals  and  using  the  normal  distribution 
for  pa(a)-*  yields: 


+  <  \ 

c  vyo)  = 


GO 

V  !'V+ 

2v  (0)  \L  i  /v  J 
C  1  ?i 
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exp(-u2)du-^  -  exp(-u2)du| 
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A  change  in  the  index  of  summation  in  the  second  sum  produces! 


-  (N.  -  N,*,) erCc? . 
1  ill  1 


(15) 


where  the  complementary  error  function  has  its  standard 
definition  (ref.  20)  : 


erfc (x) 


co 

2  P  2 

5  -j-  (  exp(-u  )  du 

x 


(IS) 


Now  note  that  -(N^-N^+j)  is  plus  one  when  is  a  positive 
maximum  or  negative  minimum  and  minus  one  if  is  a  positive 
minimum  or  negative  maximum.  This  leads  to 


The  argument  to  be  presented  is  obviously  valid  for  other  im¬ 
pulse  strength  distributions;  however,  use  of  the  normal  dis¬ 
tribution  for  pa(a)  provides  a  good  statistical  model  for  many 
environmental  noises  and  has  been  extensively  used  in 
empirical  studies  including  those  here  reported. 
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c+(yQ) 


v  i  v 

2vc(0)  l  L 

pos  max 
neg  min 


erfc!  .  - 
i 


erfcsj.  (17) 


neg  max 
pos  min 


It  is  customary  in  data- reduction  to  normalize  the  level  yQ 
with  respect  to  the  rms  of  the  response  process,  oy.  This 
rms  level  follows  by  an  application  of  Campbell's  theorem 
(ref.  6)  s 


2  5  P  2 

oy  -  vo4  j  h  ( t)  dt. 


Thus  the  normalized  crest  distribution  becomes: 


c+(p)  = 


JTi’oU  I  erfc[(|J'hJ(t)dt)l/2 


pos  max 
neg  min 
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erfc^(^j  h2(t)dt)1,/2 


neg  max 
pos  min 


Similar  computations  for  the  positive  maxima  and  rises 
(ref.  3)  yield  their  respective  exceedance  probabilities: 


M+(h)  = 


v  v  ,  r  ,  v  p  2  1/2  _tj _ -|# 

2vh(0)  l  ®lfCL  (2J  h  (t)dt)  |M.|J' 
pos  max  o 

neg  min 

oo 

v  v  ,  r  ,  „  r  ,2.  ,  .  ,  1/2  t)  -i 

*rfcL(2v  J  h  <c,dt’ 


Obviously  the  method  outlined  above  is  restricted  to 
the  class  of  rather  elementary  systems  whose  impulse  response 
functions  for  at  least  the  extrema  of  h(T)  and  jp°h2(t)dt  \  arc 
sufficiently  known  or  can  be  empirically  determined.  ThiS 
limitation  is  overcome  for  certain  engineering  predictions 
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in  Section  IV  where  the  results  predicted  by  Eq.  (9),  or 
those  obtained  by  specializing  Eqs .  (19) -(21)  to  the 
unimodal  system  of  Eq.  (5),  are  used  to  (approximately) 
predict  response  statistics  for  more  general  systems. 

Before  concluding  this  discussion,  the  above  mentioned 
results  of  Eqs.  (19)-(21)  for  the  1  DOF  system  will  be 
presented.  For  the  unimodal  impulse  response  of  Eq.  (5) 
the  maxima  and  minima  are  easily  calculated  yielding  for 
the  1  DOF  peak  statistics: 
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and 
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with  j"x"J  the  greatest  integer  function  x-1  <  fxl  <  x  and  y 
is  defined  in  Eq .  (10) .  Although  it  is  not  immediately  obvious 
from  the  form  of  these  expressions,  for  large  0  the  statistics 
depend  only  on  the  impuls ive- noise-dens ity  parameter  y  while 
exhibiting  negligible  dependence  on  Q,  v,  or  u>Q  individually. 
This  agrees  with  the  result  of  Eq.  (9)  derived  on  the  basis 
of  a  high  system  Q.  The  physical  significance  of  the  system 
parameter  o^q/Q  will  be  discussed  in  Section  IV. 


2.  GENERA LIZED-t  IMPULSIVE  NOISE  MODEL 


The  basic  hypotheses  of  Hall's  "general ized- t H  impulse 
noise  model  in  a  formal  way  facilitate  the  calculation  of  both 
crest  and  maximum  statistics.  Although  the  physical  impli¬ 
cations  of  the  final  results  caution  against  using  Hall's 
model  to  characterize  system  wear  data  one  can  nevertheless 
appreciate  its  analytical  tractability .  For  the  present 
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discussion,  the  most  significant  assumption  is  that  the 
modulating  process  m(t)  be  slowly- varying  when  compared 
to  n(t)  As  will  be  shown,  this  hypothesis  is  most  crucial. 

Consider,  then,  the  process  characterized  in  Eq.  (2)  . 
Differentiating  Eq.  (2)  twice  yields 


y(t)  *  m(t)n(t)  +  ra(t)n(t) 


(24) 


and 


y(t)  -  m(t)n(t)  +  2if\(t)n(t)  +  m(t)n(t)  .  (25) 


If  m(t)  varies  slowly  with  respect  to  n(t),  then  only  the 
first  term  of  the  right-hand  members  of  Eqs.  (24)  and  (25) 
need  be  retained  (ref.  7) .  One  can  use  this  simplification 
to  evaluate  Eq.  (7).  The  density  function 
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(26) 


follows  from  a  change  of  variables  in  Pn,n,m^"  ••  •)  to  obtain 

Pan. tan. ml . )  *  Py,v,m( . >  *nd  then' integrating  over  all 

values  of  m.  Since  n(t)  is  Gaussian  and  is  independent  of 
n(t).  one  obtains 


P  (m)  2 
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p  •  (n,n,m)  =  — - exp4  - - r 

*n, n,m  2ro  o.  c l  ,2 
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(27) 


where  (assuming  n(t)  is  zero-mean) 

o  2  *  R  (0) 
n  n 


o  • 
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d  R  (x) 
n 

dT2 


T-0 


(28) 


^Otherwise  the  "peak"  statistics  get  rather  complicated.  After 
all,  any  process  can  be  represented  as  a  Gaussian  n(t)  times 
some  "finagle  factor"  process  m(t) . 
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Combining  Eqs .  (7),  (26),  and  (27)  and  performing  the 
integration  over  the  variable  ^  leads  to  the  normalized 
crest  statistics  , 


c+(y  )  ■»  dm  p  (m)e  n  .  (29) 

o  J  m 

-00 

Thus,  the  crest  statistics  for  Hall's  model  are  obtained  in 
quadrature.  It  is  also  evident  that  appropriate  choice  of 
may  lead  to  closed-form  evaluation  of  Eq.  (29)  . 

An  approach  similar  to  that  used  for  the  derivation  of 
Eq.  (29)  also  can  be  employed  to  derive  an  expression  for 
the  maximum  statistics.  Using  the  "slowly  varying" 
assumption  on  m(t)  again  simplifies  Eqs.  (24)  and  (25), 
allowing  the  density  function  Py , y, y (•<•». )  to  be  obtained  from 
Pn,  A,  K,  i-n  a  fashion  similar  to  that  used  for  the 

crest  statistics.  The  general  result  for  the  maximum  sta¬ 
tistics  is  quite  formidable  unless  a  narrowband  assumption 
is  made  for  n(t) .  For  this  case  one  finds  that  the  maxima 
exceedance  distribution  M+(yc)  is  identical  to  c+(yD)  and 
thus  is  given  by  Eq.  (29) .  if  n(t)  is  not  narrowband,  then 
wMyo)  i-8  more  complicated.  Nevertheless,  it  is  in  quadrature 
and  involves  explicitly  only  the  exponential  and  the  error 
function.  The  density  function  pjn(*)  appears  within  the 
integrand  as  a  multiplier,  just  as  in  Eq.  (29) .  Since  the 
maximum  statistics  of  a  multimodal  linear  system  are  not 
examined  experimentally  in  this  report,  the  main  point  of 
this  discussion  is  primarily  to  illustrate  the  analytical 
tractability  of  Hall's  model. 

The  results  of  Eq.  (29)  and  the  discussion  above  of  the 
response  maximum  statistics  may  be  obtained  in  another 
"intuitive"  fashion.  One  might  envision  the  product 
representation  of  Eq .  (2)  as  a  "scaling"  of  the  Gaussian 
variate  n(t)  by  an  independent  slowly-varying  "scale"  factor 
m(t).  Thus,  the  results  for  both  the  crest  and  maximum 
statistics  of  y ( t)  can  be  obtained  by  averaging  with  weight 
Pmfm)  both  the  crest  and  maximum  statistics  of  n(t.)  over  the 
values  of  m(t),  which  may  be  regarded  formally  as  variations 
in  the  rms  level  of  n(t) .  Such  a  procedure  is  used  in  the 
studies  of  atmospheric  gust  loading  (refs.  21,  22) .  In 
applications,  however,  this  approach  must  be  used  with 
caution,  since  its  validity  is  ensured  only  when  m(t)  is 
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indeed  slowly  varying  compared  with  n(t).  Otherwise,  mort« 
information  about  the  simultaneous  statistical  behavior  of 
m(t)  and  its  two  derivatives  must  be  known. 

Whichever  analysis  is  used,  one  observes  that  numerical 
results  for  the  "peak"  distributions  depend  critically  upon 
part  icular  choices  for  the  density  function  pm(m) .  A 
particular  choice  of  Pm(m)  was  made  by  Hall  for  modeling 
of  atmospheric  radio  noise  data  (ref.  23);  for  such  noise, 
he  argued  that  a  chi-squared  distribution  for  the  random 
variable  b  £  1/m  should  be  chosen  (Eq.  (4)) .  If  so,  Eq.  (29) 
has  a  closed- form  expression,  namely  (using  Hall's  notation) 


c+(yQ) 


-y  e^2 

rl  ♦  <^)  ] 2 


(30) 


where 


9  =  M  +  1  . 


(31) 

(32) 


As  mentioned  earlier  numerical  choices  of  values  for  the 
parameters  (0, yH)  in  a  particular  study  must  await  fitting 
of  Eq.  (30)  to  observed  data. 


Rather  than  assuming  an  expression  for  p^tm)  ,  one  might 
infer  its  form  from  observed  data.  This  can  be  accomplished 
by  rewriting  Eq.  (29)  to  obtain 


c+(yD) 


CD 

f* 

f  (  t)  e 
J 
o 


(33) 


where 


f St)  -  Pni„-1/2, 


(34) 
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Now  £ q .  (33)  can  be  recognized  as  a  Laplace  transform 
( re  f  .  24  )  : 


c+(y  )  -X  !"f  (  t)  (s) 

O  L  w 


(35) 


Therefore,  if  empirical  crest  data  are  fitted  with  some 
analytical  expression,  then  f(t)  is  determined  by  the  inverse 
Laplace  transform  of  the  expression5  for  c+(yQ) ,  In  principle 
this  would  lead  via  Eq .  (34)  to  the  determination  of  p^m)  . 
Such  a  technique  would  rely  heavily  on  one’s  ingenuity  in 
finding  an  analytical  expression  c+(yQ),  yielding,  nevertheles 
an  appropriate  form  for  Pm(-)  •  In  turn,  knowledge  of  Pjn(') 
could  be  used  to  predict  other  response  statistics,  e.g., 
maxima  distributions  or  first-order  probability  densities 
which  could  be  used  to  check  the  "slowly-varying"  aspect  of 
the  model  empirically.  It  also  should  be  noted  at  this 
point,  that  for  large  arguments,  Eq .  (30)  behaves  as  an 
inverse  power.  Thus,  if  Eq .  (30)  is  to  be  used  as  a  "peak*' 
exceedance  distribution  to  calculate  partial  damage  moments, 
one  can  expect  the  higher-order  moments  to  diverge. 


^This  approach  does  not  work  too  well  with  the  Wt ibull  fits 
to  "peak"  data  reported  earlier  (ref.  5),  since  this  dis¬ 
tribution  is  seldom  transformable  in  closed  form. 


IV 


ENGINEERING  APPROXIMATION:  GENERALIZATION  OF  y 


The  preceding  sections  have  dealt  with  properties  of 
formal  analytical  solutions  for  the  response  crest  and  maximum 
statistics.  Unfortunately,  these  results  are  fairly  complicated 
for  nongaussian  noise  models.  Also,  for  general  systems 
analysis  the  reliability  engineer  may  not  know  or  be  able 
to  determine  even  such  basic  information  as  the  system 
impulse  response  function  h(t) .  The  work  of  Section  III-l 
reveals  that  only  knowledge  of  the  values  of  the  "peaks" 
of  the  system  impulse  response  is  absolutely  necessary  for 
predicting  "peak"  statistics  in  the  low  impulse  limit. 

However,  even  this  information  may  be  inaccessible  at  least 
before  the  system  has  been  built.  One  might  wonder,  then, 
if  some  simple,  statistically  reliable  characterization  of 
the  linear  system  under  study  would  not  suffice  for 
"engineering"  prediction  of  the  response  "peak"  distributions. 

In  this  section  two  approaches  (time-domain  and  frequency- 
domain),  both  based  on  considerably  less  system  information 
than  such  full  knowledge  about  h(t),  are  discussed  and  applied 
to  predicting  response  crest  statistics.  Basically,  both 
methods  consist  of  modeling  a  complex  multimodal  system  by 
an  appropriate  unimodal  system  chosen  in  such  a  way  that  the 
level  crossings  of  this  “equivalent"  system  response  approxi¬ 
mate  those  of  the  complex  system. 

The  basis  for  this  system  modeling  lies  in  the  simple 
dependence  cf  the  unimodal  statistics  on  certain  system 
parameters.  As  cited  earlier,  for  a  1  DOF  system,  y  =  QvAuo 
completely  determines  the  crest  statistics.  The  system 
characterization  reflected  by  y  appears  only  as  the  ratio 
Q/<dq.  Thus,  in  order  to  model  a  complex  system  by  a  unimodal 
system  in  a  manner  appropriate  to  predicting  level  crossings 
it  is  first  necessary  to  determine  a  more  general  system 
parameter,  analogous  to  Q/“0-  for  the  unimodal  system  which 
in  turn  can  be  applied  to  arbitrary  systems  with  an  equivalent 
positive  slope  zero-crossing  rate  cOy/27:.  These  two  parameters, 
then,  determine  the  crest  statistics  for  the  "equivalent" 
unimodal  system.  It  is  logical  to  ask  if  the  crest  statistics 
predicted  in  Section  III-l  for  a  unimodal  system  can  be  used 
to  obtain  an  approximation  to  the  complex  system  crest 
statistics.  This  possibility  will  now  be  explored. 
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1. 


TIME-DOMAIN  APPROACH 


In  addition  to  lending  some  weight  to  our  hypothesis  of 
the  validity  of  a  generalized  definition  of  y  as  outlined 
above,  Middleton's  work  (ref.  8)  on  Poisson  impulsive  noise 
suggests  a  method  of  making  a  meaningful  extension  of  the 
unimodal  parameter  Q/^>0  to  complex  systems.  His  results 
predict  that,  for  low  impulse  frequencies,  the  nth  order 
response  probability  densities  (and  thus  the  frequency  of 
level  crossings  via  Eq.  (7))  depend  strongly  on  the  product 
of  the  impulse  frequency  v  and  the  duration  of  the  system's 
impulse  response,  t.  This  agrees  exactly  with  our  unimodal 
crest  results,  since  0/,JJo  is  directly  proportional  to  the 
time  required  for  the  exponential  impulse  response  envelope 
to  decay  to  a  given  value  and  QvAu0  determines  the 
system  crest  statistics.  The  problem  now  is  to  select  a 
convenient,  physically  meaningful  measure  of  impulse  response 
duration  for  a  more  complex  system. 


The  most  convenient  definition  of  response  duration 
may  well  depend  on  the  particular  system  under  study. 
Initially  we  restrict  our  attention  to  linear  multimodal 
systems  having  a  fairly  high  number  of  modes  since  the 
empirical  results  of  Section  V  deal  with  such  systems. 

For  the  statistical  analysis  of  room  acoustics 
acousticians  (ref.  25)  have  found  it  convenient  to  describe 
the  room  (system)  in  terms  of  its  reverberation  curve  (the 
time  history  of  sound  pressure  level  following  cessation  of 
white  noise  excitation).  In  particular,  the  system  rever¬ 
beration  time,  the  60  dB  decay  time  based  on  an  exponential 
fit  of  the  initial  portion  of  the  reverberation  curve,  pro¬ 
vides  a  convenient  measure  of  the  duration  of  the  impulse 
response  for  multimodal  systems  in  general.  In  order  to  have 
application  here  this  "reverberation  time",  'tfc0dB'  must  be 
converted  to  an  e-^/^  decay  time  t  in  conformity  with  the 
definition  of  QAdq  as  a  characteristic  decay  time  for 
the  1  DOF  system.  Thus  the  e" 

T60dB  as 


1/2 


decay  time  t  is  related  to 


60  dB 

6  loa  10 
e 


(36) 


Using  i  as  defined  above  as  the  generalization  of  Q/u>0  (i.e., 
Y  =  vt)  yields  the  crest  statistics  of  the  "equivalent" 
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unimodal  system  via  Section  III-l.  Methods  of  approximating 
the  complex  system  response  on  the  basis  of  th's  unimodal 
system  response  will  be  discussed  in  the  section. 

To  date  no  rigorous  analytical  justification  has  been 
found  for  the  choice  of  the  time  constant  t  obtained  from  the  \ 

system's  reverberation  curve  as  a  meaningful  measure  of  j 

response  duration.  Empirical  verification  aside,  the  1 

reverberation  time,  being  based  on  the  rate  of  system  energy 

decay  in  a  manner  analogous  to  QAt>0  for  the  1  DOF  system,  j 

seems  to  be  intuitively  reasonable  as  a  choice  for  generalizing 
Y  in  terms  of  a  time-domain  parameter.  However,  its  ultimate 
justification  must  lie  in  empirical  verification. 

2.  FREQUENCY  DOMAIN  APPROACH 

Very  often,  due  to  the  nature  of  the  system,  it  is 
preferable  to  work  with  frequency-domain  parameters  rather 
than  to  use  the  time-domain  approach  outlined  above.  In 
any  case,  the  possibility  of  a  generalized  interpretation 
of  y  carries  over.  The  relationship  existing  between  time- 
domain  and  frequency-domain  concepts  merely  converts  the 
impulse  response  duration  information  into  a  certain  measure 
of  system  bandwidth  B.  As  was  the  case  with  response 
duration,  t ,  system  bandwidth  is  a  rather  vague  term.  There 
exist  a  variety  of  definitions  for  bandwidth  in  the  litera¬ 
ture  (ref.  26)  each  defined  in  a  manner  useful  for  the 
particular  problem  area  involved.  For  example  the  "equivalent 
noise  bandwidth"  defined  by: 

oo 

c  2 

.  |  h  (  ja>)  |  cU 

u 
o 

B  *  - 3 - 

n  |H(ju))| 

max 

with  H ( juj)  the  fourier  transform  of  h(r),  is  useful  in 
describing  the  noise  power  transmitted  by  a  systc..)  under 
white  noise  excitation,  while  the  "3  dB  bandwidth",  defined 
as  the  difference  between  the  two  frequencies  and  c vp  at 
which  |  H  ( juj)  |  ^  is  3  dB  below  its  peak  value,  is  convenient 
for  describing  relatively  flat  bandpass  systems.  Most 
definitions  for  B  that  are  applicable  to  the  unimodal  system 
yield  a  bandwidth 


I 

i 

i 

(37)  ! 
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(38) 


For  more  complex  systems,  however,  the  various  bandwidth 
expressions  differ,  each  depending  strongly  on  which 
particular  aspect  of  the  system  transfer  characteristic  is 
emphas ized . 


Our  immediate  objective  is  to  find  a  definition  of 
bandwidth  that,  is  a  measure  of  the  inverse  of  the  aforementioned 
response  duration  for  the  highly  peaked  multimodal  response 
spectra  of  interest  in  structural  vibrations.  Let  us  now 
return  to  the  system  time-domain  response.  For  the  type  of 
system  under  consideration,  response  decay  is  governed  by 
the  rate  of  system  energy  dissipation.  Thus,  a  bandwidth 
expression  based  on  such  energy  considerations  is  sought. 

To  relate  this  to  the  1  DOF  system,  consider  the  fundamental 
definition  of  Q  (ref.  27): 


2<o  T 


where  T  is  the  system  time-average  kinetic  energy  and  D  is 
its  time-average  rate  of  power  dissipation.  Combining 
Eq .  (39)  with  Eq .  (38)  gives  the  unimodal  bandwidth  expression 
in  terms  of  the  systrm  energies  as: 


However,  the  quantities  D  and  T  are  well  defined  even  for 
arbitrary  systems;  hence,  Eq .  (40)  will  be  used  as  our 
general  bandwidth  definition,  being  a  measure  of  the  rate 
of  decay  ot  system  energy. 


To  obtain  useful  results  from  this  expression  let  us 
specialize  our  attention  to  multimodal  systems  having  highly 
peaked,  well  defined  modes  such  as  the  circular  metal  plate. 

If  the  modes  are  of  sufficiently  high  Q,  they  can  be  assumed 
to  be  uncoupled,  justifying  writing  the  total  system  kinetic 
energy  as  the  sum  of  the  individual  modal  energies  T^  (ref.  28) . 
Similarly  the  total  rate  of  energy  dissipation  is  the  sum  of 
the  individual  rates  Dy.  Now,  for  each  mode,  Eq .  (40)  can 

be  applied,  again  assuming  high  modal  Q's,  to  express  Dj 
in  terms  of  T^,  x  =  2BiTi  with  Bx  the  modal  3  dB  bandwidths . 
Thus  ,gq .  (40)  for  the  multimodal  sys  ten.  becomes 
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(41) 


V  B  .  T  . 
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That  is,  the  system  bandwidth  is  just  the  average  of  the 
individual  modal  3  dB  bandwidths  weighted  with  respect  to 
the  modal  kinetic  energies.  A  more  general  interpretation 
of  y,  based  on  Eq .  (40)  becomes 


A  certain  amount  of  care  must  be  exercized  in  applying 
Eq.  (41)  to  distributed  systems  in  general.  In  particular, 
when  considering  response  statistics  concern  is  solely  with 
a  particular  system  input-output  characteristics,  i.e.,  the 
response  at  a  single  fixed  point,  and  not  with  over-all 
system  behavior.  From  a  reliability  viewpoint  one  is 
concerned  only  with  response  t  a  few  selected  points  of 
greatest  system  vulnerability.  Hence,  each  modal  energy 
must  not  be  interpreted  as  an  over-all  modal  energy  summed 
spatially  over  the  eigenmode  shape,  but  rather  as  the 
"localized’'  modal  energy  at  the  response  point  in  question, 
i.e.,  the  kinetic  energy,  due  to  the  ith  mode,  of  an  infinitesimal 
mass  element  at  the  response  point  as  predicted  by  the  actual 
system  transfer  characteristic.  If  the  Tj/s  are  interpreted 
as  modal  energies  averaged  over  the  plate,  Eq .  (41)  yields 
an  "average  system  bandwidth"  that  is  independent  of  the 
fixed  response  point  and  has  no  significance  in  so  far  as 
the  actual  system  input-output  frequency  characteristics 
for  the  localized  response  variable  under  analysis  are 
concerned . 

Even  with  the  above  two  methods  for  selecting  a  unimodal 
system  "equivalent"  to  the  given  complex  system  (or  more 
precisely  the  equivalent  system  v)  available,  the  problem 
of  determining  a  relationship  between  the  actual  system 
crest  statistics  and  those  of  the  unimodal  system  remains, 

Since  the  goal  here  is  to  obtain  rather  simplified  engineering 
approximations,  the  entire  development  of  this  section  lacks 
a  rigorous  analytical  basis,  and  any  justification  for  the 
approximations  must  lie  with  experiments.  Similarly,  the 
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relationship  between  the  statistics  of  the  equivalent 
unimodal  system  and  the  modeled  system  is  based  partially 
on  qualitative  consideration  of  the  two  system  impulse 
responses  and  ultimately  on  empirical  results.  This 
rela t ions'n ip  will  be  investigated  next. 

The  moat  convenient  relationship  between  the  unimodal 
and  complex  multi-modal  system  statistics,  if  it  yields  good 
approximations,  is  a  direct  equality.  While  a  considerable 
amount  of  empirical  data  on  various  multimodal  systems 
(cf.  Section  V)  indicates  that  the  functional  form  of  the 
unimodal  crest  distribution  agrees  with  measured  multimodal 
results,  in  most  cases  the  two  differ  by  a  multiplicative 
factor  of  the  order  of  unity.  This  discrepancy  is  explained 
qualitatively  by  consideration  of  the  impulse  responses  of 
the  two  systems  involved.  The  unimodal  response,  being 
sinusoidal  in  nature,  has  a  zero-crossing  preceding  each 
non-zero  level  crossing.  For  a  multimodal  system,  on  the 
other  hand,  in  most  cases  the  modes  add  in  such  a  manner 
that  there  may  be  several  extrema  of  the  response  (producing 
level  o;  ossings  at  levels  i|  >  0)  between  successive  zero- 
crcseings.  Thus  a  unimodal  system  having  the  same  frequency 
of  level  crossings  at  large  response  levels  as  the  multimodal 
system  cannot  have  the  correct  level  crossing  rate  for  small 
response  levels  (say  v  ~  0) .  Hence  while  the  unimodal  system 
selected  as  indicated  in  Sections  IV-1  or  IV-2  correctly 
predicts  the  frequency  of  level  crossings  for  the  modeled 
system  for  t)  >  0.5,  it  cannot  predict  the  normalized  crest 
distribution.  This  is  because  c+(g)  is  obtained  by  normalizing 
the  predicted  1  ’el  crossing  frequency  with  respect  to  the 
equ  ivalent  unimodal  resonant  frequency  (zero-crossing  rate) 
uj^/2-n- ,  On  the  other  hand  the  complex  system  normalized 
crest  statistics  are  obtained  by  normalizing  the  level 
crossing  frequency  with  respect  to  the  true  system  zero¬ 
crossing  rate,  which  is  slightly  smaller  than  the  zero-crossing 
rate  for  the  unimodal  system,  o>0/2'n,  for  the  reasons  stated 
above . 

Thus,  the  crest  predictions  differ  from  the  true 
statistics  by  a  multiplicative  factor,  the  ratio  of  the  two  zero- 
crossing  rates  which  predicts  a  larger  "peak"  exceedance 
probability  at  higher  levels  t.hon  for  the  unimodal  case.  Hence, 
it  appears  that  in  normalizing  one  should  measure  or  estimate 
the  ac tua 1  frequency  of  zero-crossings  for  the  complex  system 
rather  than  use  the  zero-crossings  corresponding  to  the 
equivalent  unimodal  frequency.  Examination  of  the  derivations 
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in  Section  Ill-l  reveals  that  the  level  crossing  frequency 
for  the  unimodal  system,  Vy+(ri),  is  given  by 

u> 

VY+<'')  =  2?  Cy+(T1)*  (43) 

Here  cy+(ti)  is  the  normalized  crest  distribution  for  the 
unimodal  system  with  impulse  density  parameter  y.  Using 
Eq.  (43)  as  the  predicted  level  crossing  frequency  for  the 
complex  system  and  normalizing  with  respect  to  the  actual 
(ox  estimated)  complex  system  zero-c tossing  rate,  v+(0), 
gives  the  approximate  crest  distribution,  c+(h),  for  the 
complex  system: 

CD 

c+  (  r\)  = - 7 -  c  +(Ti)  (44) 

2ttv  (0)  Y 

where  y  is  obtained  by  the  methods  given  at  the  beginning  of 
this  section  and  the  equivalent  u>0  is  yet  to  be  determined. 

In  a  development  similar  to  that  of  Eqs .  (39) -(42) 
for  multimodal  systems,  Ungar  (ref.  28)  suggests  for  this 
equivalent  resonant  frequency,  u>0,  the  center  frequency  of 
a  band  encompassing  all  modes  of  the  system  that  give  signi¬ 
ficant  contributions  to  the  bandwidth  expr  ssion  Eq4  (41) . 
Thus  the  band  is  chosen  to  contain  essentially  all  f of  the 
system  energy.  For  the  highly  peaked  spectra  studied  here 
this  band  is  well  defined;  for  other  systems  this  might  not 
be  the  case.  A  possible  alternative  expression  for  is 
discussed  in  Section  VI.  However,  since  the  center  frequency 
expression  above  yields  good  agreement  with  empirical  statis¬ 
tics  for  multimodal  systems  it  is  used  for  all  the  results 
presented  in  Section  V. 

The  two  methods  outlined  here  for  selecting  the  equi¬ 
valent  system  y  (and  to0)  are  those  appropriate  to  predicting 
the  lightly-damped  circular  plate  response  statistics,  i.e., 
methods  relying  on  the  highly-peaked,  multimodal  nature  of 
the  system  transfer  characteristic.  For  other  types  of 
systems,  while  the  exact  expressions  derived  here  may  not  be 
applicable,  the  concept  of  selecting  an  equivalent  y  and  u> 
may  still  be  applied,  perhaps  with  scmt  other  definition  of 
/stem  bandwidth  or  response  duration  more  appropriate  to  the 
system  under  consideration. 
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V 


COMPARISONS  WITH  EMPIRICAL  CREST  DISTRIBUTIONS 


Several  methods  have  been  developed  in  the  preceding 
sections  which  should  enable  one  to  predict  with  reasonable 
accuracy  the  response  crest  statistics.  In  this  section 
the  application  of  both  the  "exact"  theories  and  the 
"equivalent  unimodal  system"  methods  will  be  discussed. 

In  particular,  empirical  crest  data  obtained  from  impulsive 
source  excitation  of  a  clamped  circular  plate,  along  with 
data  obtained  from  unimodal  and  bimodal  systems,  will  be 
compared  with  the  predictions  mentioned  above. 

1.  EXPERIMENTAL  EQUIPMENT 

The  analytical  impulsive  noise  model  of  Eq .  (1)  is 
obviously  not  realizable  by  a  physical  noise  generator.  As 
reported  in  previous  work  (ref.  10)  a  generator  approximating 
this  idealized  model  has  been  developed.  This  generator 
produces  a  train  of  approximately  Poisson-distributed,  fixed- 
duration,  rectangular  pulses  whose  amplitude  distribution 
pa(a)  is  fixed  by  modulating  the  pulse  train  with  a  random 
noise  having  the  desired  probability  density.  For  all  data 
presented  in  this  section  pa(a)  was  the  normal  density, 
with  the  pulse  width  T  fixed  such  that  if  is  much  greater 
than  the  greatest  response  frequency  of  the  system  to  be 
excited,  the  generator  produces  a  good  approximation  to  the 
ideal  noise  model  of  Eq.  (1)  . 

Data  obtained  using  this  generator  will  be  presented  for 
three  types  of  systems:  unimodal  (.1  DOF)  systems,  bimodal 
systems  and  the  (multimodal)  circular  plate.  The  unimodal 
systems  were  obtained  by  synthesizing  the  equation: 

y  +  — -  y  +  o>o  y  =  f(t)  (45) 


with  analog  computer  modules  (ref.  11) .  The  forcing  noise 
is  denoted  by  f(t).  This  simple  system  is  included  only 
to  verify  Eqs .  (9)  and  (22) .  The  second  type  of  system 
was  a  bimodal  system  consisting  of  a  pair  of  elastically 
coupled  oscillators,  again  synthesized  with  analog  computer 
modules  (ref.  3).  Specifically,  the  equations  synthesized 
were ; 
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f  (t) 


mnyi  +  Cliyi  +  kUyi  +  k12y2  - 

m22y 2  +  C22y2  +  k22y2  +  k2iyl  =  gf(t) 


(46) 


wtch  yj  the  observed  response  variable  and  g,  n\jk,  Cj^  and  k-k 
positive  constants.  In  general  these  equations  produce  a 
system  transfer  characteristic  exhibiting  two  distinct 
resonant  peaks.  Various  relative  mode  strengths  were  obtained 
by  varying  the  Q's  of  the  two  oscillators.  Peak  statistics 
obtained  on  these  two  types  of  systems  have  been  reported 
previously  (refs.  3,  5) .  The  major  emphasis  here  will  be 
on  a  complex  multi-modal  system,  the  clamped  circular  plate. 
Details  of  the  plate  characteristics  are  included  in  the 
Appendix. 

2.  "EXACT  THEORIES" 

a.  Poisson  Impulsive  Noise  -  Low-v  Case 

The  analytical  methods  ot  Section  III-l  can  be  applied 
to  predicting  the  response  crest  distribution  of  a  unimodal 
system.  For  these  purposes,  Eq.  (9)  and  (22)  apply.  Direct 
comparisons  reveal  that  the  more  general  result  of  Eq.  (22) 
agrees  in  detail  with  the  earlier  published  results. 

Figure  2  presents  response  crest  statistics  for  nine  values 
of  the  impulse  noise  density  parameter  y;  the  data  can  then 
be  interpreted  as  either  experimental  results  or  the 
predictions  of  Eq.  (9)  and  Eq .  (22);  all  three  being  the 
same  for  a  1  DOF  system  (only,  of  course,  for  y  <  0.5). 

The  gaussian  limit,  y  ■*  oo ,  is  shown  for  comparison  purposes. 

One  is  encouraged,  then,  to  test  Eq.  (19),  from  which 
Fq ,  (22)  is  derived,  for  more  complex  systems. 

The  bimodal  system  provides  a  convenient  system  for 
this  te9t  (ref.  3)  .  However,  even  for  this  system  the  form 
of  the  impulse  response  is  not  readily  adapted  to  the  analysis. 
Thus  to  evaluate  Eq .  (19),  the  maxima  and  minima  of  h(r) 
were  measured  empirically.  With  the  individual  modal  Q's 
set  at  various  values  in  the  range  4  <  Q  <  20  it  is  necessary 
to  retain  from  ten- to- twenty  terms  in  the  summation  of  Eq.  (19) . 
Data  for  a  typical  system  is  compared  to  the  prediction  of 
Eq.  (19)  in  Fig.  3  for  four  values  of  v.  Some  details  of 
the  transfer  characteristic  for  this  system  are  illustrated 
in  Fig.  4.  The  same  general  agreement  was  observed  for  other 
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bimodal  system  transfer  characteristics  (ref.  3),  namely 
good  agreement  for  the  lower  impulse  frequencies  with  some 
deviation  for  the  highest  frequency.  This  discrepancy  will 
be  discussed  further  in  Section  VI.  It  should  be  noted 
however,  even  for  the  highest  impulse  frequency,  that  the 
prediction  of  Eq.  (19)  is  still  far  superior  to  assuming 
Gaussian  statistics,  which  is  currently  the  only  alternative 
for  large  impulse  frequencies. 

It  becomes  increasingly  difficult  to  apply  the  result 
of  Eq.  (19)  to  systems  more  complex  than  the  bimodal  case 
examined  above.  Consider,  for  example,  the  response  crest 
statistics  of  a  clamped  circular  plate  forced  with  impulsive 
noise.  Here  the  impulse  response  for  such  a  system  is 
extremely  complicated,  and  one  has  no  detailed  knowledge 
of  h(t)  .  in  this  case,  Eq.  (19)  would  be  inapplicable  and 
other  methods  for  examining  the  response  crest  statistics 
must  be  found. 

b.  Generalized-t  Impulsive  Noise  Model 

The  crest  statistics  of  a  noise  process  charac ter ized 
by  Hall's  model  (ref.  7)  (Section  III-2)  can  be  readily 
computed.  In  spite  of  the  shortcomings  of  his  model 
(primarily,  the  inability  to  determine  a  priori  the 
parameters  vH  and  0)  it  may  nevertheless  serve  as  a  starting 
point  for  determining  an  analytical  form  for  the  response 
crest  statistics  of  a  particular  multimodal  system. 

To  this  end,  then,  yH  and  0  can  be  determined  by  fitting 
Eq .  (30)  to  some  typical  response  data.  This  task  is  better 
accomplished  after  first  noting  that  the  ultimate  application 
ia  the  calculation  of  the  system  damage  moments.  As  a  result 
it  will  probably  be  most  important  to  fit  the  data  for  the 
larger  response  levels  where  the  contribution  to  these  moments 
is  most  important.  Figure  5  indicates  the  results  of  such  an 
operation  for  three  typical  cases  of  multimodal  (i.e., 
clamped  circular  plate)  response  crest  statistics.  The 
various  system  transfer  characteristics  of  these  three 
typical  cases  are  shown  in  Fig.  6. 

The  agreement  between  Eq .  (30)  and  multimodal  experimental 
data  is  quite  remarkable.  Although  some  discrepancies  do 
exist  at  the  lower  response  levels,  they  make  little  dif¬ 
ference  insofar  as  damage  calculations  are  concerned.  And, 
in  fact,  since  the  fitted  curve  seems  to  lie  above  experimental 
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data  for  these  lower  response  levels,  the  analytical  moments 
should  provide  a  conservative  estimate  of  the  damage.  By 
and  large  such  a  "fitting"  procedure  parallels  earlier 
attempts  to  describe  response  data  with  Weibull  distributions 
(ref.  5).  However,  unlike  the  case  of  the  Weibull  distri¬ 
bution,  the  evaluation  of  yH  and  0  here  will  critically 
delimit  the  convergence  of  even  lower-order  damage  moments 
which  "physically"  are  known  to  be  finite.  Thus,  more 
appropriate  methods  of  characterizing  the  response  crest 
statistics  of  multimodal  systems  are  sought. 

3.  EQUIVALENT  UNIMODAL  APPROXIMATIONS 

Even  though  the  exact  results  of  Section  III  yield  very 
good  predictions  of  general  system  peak  statistics,  in  many 
systems  they  do  not  apply.  Even  for  a  bimodal  system  their 
applicability  generally  relies  on  an  empirical  determination 
of  h(t) .  For  systems  such  as  the  high-Q  circular  plate, 
empirical  measurements  of  h(T)  are  impractical.  Thus,  in 
this  section,  the  approximate  methods  derived  in  Section  IV 
will  be  applied  to  the  circular  plate  and  to  the  bimodal 
systems  of  Section  V-l. 

For  this  application,  it  is  necessary  to  determine 
either  the  system  reverberation  time  t  or  the  system  bandwidth 
of  Eq.  (41)  in  order  to  apply  the  results  of  Section  IV. 

For  the  circular  plate  the  reverberation  time  based  on  the 
initial  decay  slope  was  measured  using  standard  techniques 
(ref.  29)  with  a  Bruel  and  Kjaer  level  recorder  ana  noise 
generator.  For  reasons  to  be  discussed  in  Section  VI  this 
equipment  was  not  suitable  for  measurements  of  the  reverbera¬ 
tion  time  of  the  bimodal  system.  To  date  no  acceptable  alterna¬ 
tive  method  has  been  found  for  measuring  decay  time  of  a 
bimodal  system.  Hence,  all  predictions  will  be  based  on  the 
system  bandwidth,  Eq .  (41)  .  For  the  circular  plate,  results 
will  be  presented  based  both  on  the  reverberation  time  and 
the  system  bandwidth. 

Before  presenting  the  data  let  us  consider  the  bandwidth 
prediction.  In  evaluating  the  bandwidth  of  the  circular 
plate  from  Eq.  (41)  it  is  necessary  to  determine  the  modal 
kinetic  energies  at  the  response  point.  it  is  obvious  from 
a  typical  system  transfer  characteristic  (cf.  Fig.  6)  that 
the  vibrational  energy  of  a  particular  mode  of  the  plate  is 
concentrated  in  a  narrow  frequency  band  about  that  modal 
resonance  frequency  and  thus  has  negligible  response  at  any 
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other  modal  resonant  frequency.  ThuB,  to  the  degree  of 
approximation  involved  in  the  method  of  Section  IV  each 
particular  mode  will  be  assumed  to  be  a  high  Q  resonant 
peak  similar  to  the  1  DOF  system.  From  '.his  assumption 
the  modal  kinetic  energies  become: 


T  .  =  t  m  P  .  B  .ui. 
1  4  ill 


(47) 


where  is  the  absolute  response  mean  square  at  the  resonant 
frequency  tu •  and  m  is  the  surface  mass  density  at  the 
response  point.  Thus,  the  bandwidth  of  the  system  becomes: 


B  =  - - r-  .  (48) 

P  B  ,oj  . 

-  Ill 

j 

From  this  expression  one  gains  some  insight  into  the 
relative  significance  of  the  various  response  modes.  In 
particular,  for  a  typical  system,  with  B ^  roughly  proportional 
to  the  modal  bandwidths  are  weighted  approximately 

by  the  third  power  of  modal  resonant  frequency.  Thus  a  high 
frequency  mode  whose  response  is  an  order  of  magnitude  or 
more  below  the  fundamental  mode  will  have  a  much  more 
significant  contribution  to  system  bandwidth  than  will  the 
fundamental.  Indeed,  for  the  system  of  Fig.  6-b  the  contri¬ 
bution  of  the  fundamental  mode  to  Eq .  (48)  is  insignificant 
compared  to  those  of  the  sixth  and  seventh  modes  whose 
responses  are  20  dB  lower. 


Several  multimodal  transfer  characteristics  for  the 
circular  plate  were  studied,  obtained  by  using  various  pickup 
probe  positions  and  various  amounts  of  vibration  damping 
tape.  Selected  crest  data  so  obtained  (Propagation  Research 
Laboratory  impulsive  crest  runs  numbered  117-119)  are 
presented  in  Figs.  7-9.  The  three  figures  represent  the  sys¬ 
tem  illustrated  in  Fig.  6-b  excited  at  three  different  values 
of  average  impulse  frequency  v  =  15,  5,  1.  Also  shown  in 

these  figures  are  the  crest  statistics  predicted  by  the 
reverberation  time  and  bandwidth  methods  discussed  above. 
Generally,  the  agreement  between  the  predictions  and  the 
data  is  fairly  good  in  the  range  n  >  C.5,  probably  within  a 
factor  of  two  to  three,  for  the  nongaussian  situations 
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reptesented  (y  <  0.5) .  Also,  the  agreement  obtained  via  the 
two  methods  is  roughly  the  same;  one  difference  is  that  the 
reverberation  time  predictions  favor  more  Gaussian  behavior^ 
while  the  bandwidth  predictions  are  exemplary  of  slightly 
more  nongaussian  behavior.  Thus,  the  former  seemingly 
provides  a  lower,  while  the  latter  an  upper  bound  for  reli¬ 
ability  predictions.  Both  methods  yield  good  approximate 
crest  statistics;  the  reliability  engineer  seeking  conservative 
estimates  of  system  lifetime  would  find  the  bandwidth 
prediction  more  applicable  than  the  reverberation  time 
prediction . 

The  multimodal  bandwidth  expression  Eq.  (48)  also  will 
be  used  to  predict  the  bimodal  system  crest  statistics 
although  more  drastic  approximations  are  involved  in  this 
case.  In  general,  the  Q'  s  of  the  bimodal  system  response 
peaks  are  not  sufficiently  large  to  validate  the  assumptions 
leading  to  Eq.  (48) .  Moreover,  the  response  peak  heights 
and  bandwidths  measured  from  the  system  transfer  characteris¬ 
tic  are  not  the  true  modal  parameters.  In  general  the 
approximations  involved  tend  to  weight  the  bandwidth  of  the 
high  frequency  mode  too  strongly  in  Eq.  (41) .  Since  the 
bandwidth  of  this  mode  is  genera] ly  larger  than  the  other, 

Eq.  (48)  tends  to  predict  a  larger  system  bandwidth  than 
expected.  This,  then,  corresponds  to  a  conservative 
prediction  of  the  system  crest  statistics  insofi*  as  damage 
calculations  are  concerned. 

The  approximate  crest  predictions  of  Section  IV  using 
the  bimodal  bandwidth  of  Eq .  (48)  and  u>Q  as  defined  in 
Section  IV  are  compared  with  bimodal  crest  data  in  Fig.  10. 

It  should  be  noted  that  for  all  but  the  highest  impulse 
frequency  the  deviation  between  the  data  and  the  approximate 
prediction,  is  less  than  a  factor  of  two.  Also  the  predictions 
of  Section  IV  are  slightly  too  nongaussian  using  the  bandwidth 
expression,  Eq .  (48),  as  might  be  expected. 


^I.e.,  the  "reverberation  time"  tends  to  be  too  large  (decay 
constant  too  small)  . 
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VI 


LIMITATION  OF  THEORETICAL  PREDICTION  OF  CRESTS 


Little  need  be  said  here  about  the  agreement  of  the 
"exact"  theories  of  Section  III-1  with  experiment  for  a 
unimodal  system  excited  by  impulsive  noise.  By  in  large, 
earlier  work  (ref.  3)  has  shown  Eq .  (22)  to  be  valid  for 
sufficiently  small  values  of  the  impulse  frequency  v.  It 
is  important  to  point  out,  however,  that  the  theory  in 
Section  III-l  summarized  by  Eq .  (19)  (specifically  for  a 
unimodal  system  Eq .  (22))  has  been  found  to  agree  with 
Eq .  (9)  for  values  of  the  impulse-noise-density  parameter  y. 

This  success  encouraged  a  test  of  the  more  general  "exact" 
theory  of  Eq.  (19)  by  applying  it  to  a  more  complex  system 
and  also  encouraged  attempts  to  extend  an  interpretation  of 
y  to  complex  multi-modal  systems. 

As  an  initial  step  toward  reliability  predictions  for 
complex  systems  the  bimodal  case  was  considered  (ref.  5) . 
Again,  it  has  been  shown  (ref.  3)  that  Eq.  (19)  is  a  most 
acceptable  predictor  for  the  response  peak  statistics. 

Although  it  is  most  accurate  in  the  low  impulse- frequency 
range.  Eq .  (19)  also  can  provide  estimates  of  the  response 
crest  statistics  for  larger  values  of  v.  The  real  difficulty 
lies  in  an  inability  to  explicitly  state  how  accurate  Eq .  (19) 
is  for  a  given  value  of  v  in  a  general  way.  For  the  unimodal 
case  the  parameter  y  incorporates  both  the  noise-source 
parameter  v  and  the  system  decay  time  (reciprocal  bandwidth) 
parameter  Q/o>0  into  the  prediction.  But,  for  a  multimodal 
system  "bandwidth"  is  not  a  clear-cut,  well-defined  parameter. 
It  appears,  then,  that  ’.ntil  such  a  parameter  becomes  well- 
defined  it  may  not  be  possible  to  describe  the  general  case 
in  a  straightforward  manner. 

One  should  not,  however,  at  this  point  dismiss  the 
results  of  Eq .  (19)  as  being  purely  academic.  in  the  highly 
nongaussian  situations  it  serves  for  the  present  as  an 
accurate  predictor  for  the  response  peak  statistics  of  a 
multi-modal  system.  Knowledge  of  the  values  of  the  maxima 
and  minima  of  h(i)  is  sufficient  for  the  application  of  the 
technique.  So  interpreted,  one  is  not  required  to  know 
every  detail  of  h(T).  This  technique  will  therefore  be 
most  useful,  when  one  is  confronted  with  estimating  the  time 
to  failure  of  a  particular  system. 
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Although  application  of  the  result  Eq .  (19)  is 
straightforward  for  low  impulse  frequencies  when  the  summa¬ 
tion  rapidly  converges,  empirical  data  indicate  that  Eq .  (19) 
can  be  applied  for  impulse  frequencies  slightly  larger 
than  those  assumed  in  the  derivation.  For  such  applica¬ 
tions,  however,  in  order  to  avoid  summing  over  extrema 
occurring  after  the  beginning  of  the  next  impulse  response 
it  is  necessary  to  truncate  the  summation  to  include  only 
those  extrema  occurring  before  t  =  This  truncation  proce¬ 

dure  (which  eliminates  various  low-amplitude  peaks)  is 
included  in  the  results  of  Eq .  (2.?). 

In  an  effort  to  overcome  the  need  for  detailed  knowledge 
of  h(t).  Hall's  noise  model  was  considered  in  an  attempt  to 
characterize  system  response.  For  lack  of  justification 
of  the  use  of  any  other  modulating  density  function  pm( *) , 
Hall's  goner  a  1  ized- 1  case  was  examined.  This  led  to  a 
convenient  closed  form  expression  for  the  crest  statistics, 

Eq .  (30)  .  Unfortunately,  though,  the  parameters  yH  and  9 
can  only  be  determined  by  fittiny  Eq .  (30)  to  actual  response 
crest  data.  And  then,  it  is  possible  that  infinite  damage 
predictions  may  result.  Thus  again  one  is  restrained,  and 
the  application  of  Hall’s  model  must  proceed  on  a  “case-by- 
case"  approach.  Table  I  below  illustrates  one  such  attempt. 
The  data  ranges  from  a  highly  nongaussian  case  (run  107) 
to  a  moderately  Gaussian  one  (run  37) .  The  resultant 

Table  I  values  of  Hall's  Parameters  Determined 
by  Fit  to  Experimental  Crest  Data 


ran 

6 

37 

10.0 

2.9 

11C 

'  ■> 

O  .  4L 

1  .63 

10  7 

3 .8 

1  . 12 

calculation  of  9  and  suggests  ..hat  both  parameters 
decrease  as  more  nongaussian  processes  are  encountered.  Of 
the  two,  variations  in  the  parameter  are  perhaps  of  lesser 
importance  since  they  basically  represent  a  change  of  scale 
fer  the  response  level.  Of  greater  significance  aie  the 
changes  in  9.  The  behavior  of  Eq.  (30)  for  large  yQ  is 
inversely  proportional  to  yQ  raised  to  the  (© -2)  power.  If 
Eq .  (30)  is  to  be  used  for  the  calculation  of  damage  moments 
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(see  Section  VII)  ,  then  there  appear  to  be  some  damage 
moments  which  will  no  longer  be  finite.  In  fact,  if  the 
highest-order  finite  damage  moment  is  of  order  1 ,  then 
necessarily 

■  t  <  ;  0-2  ’  (4  9) 

r-  - 

where  ■  .  is  the  greatest  integer  operation.  Thus,  if  one 
assumed  fhe  fit  of  Hall's  analytical  model  to  be  correct 
for  the  highly  nongaussian  case  of  run  10?  a  second-order 
damage  moment  does  not  exist.  And,  even  for  the  "more 
Gaussian"  run  37,  the  damage  moments  of  order  less  than 
eight  are  the  only  ones  that  are  finite'. 

This,  of  course,  is  a  physically  unrealistic  situation. 
The  point  thaL  should  be  made  is  that  although  F.q .  (30) 
can  be  adjusted  to  fit  typical  response  data  quite  nicely, 
it  is  not  indicative  of  the  true  system  response  at  the 
larger  response  levels,  say  greater  than  ten  times  root-mean- 
squaie.  One  could  apply  some  truncation  procedure  to  Eq .  (30) 
such  an  operation  would  need  to  be  systemized,  however,  after 
more  careful  study  of  each  independent  case. 

It  appears,  then,  that  insofar  as  damage  calculations 
are  coiverned,  Hall's  "general  ized- 1"  model  is  inadequate  to 
describe  the  reliability  of  linear  systems  forced  by  impulsive 
excitation.  Future  work  may  suggest  some  other  modulating 
density  function  •)  which  will  lead  to  more  acceptable 
analytic  response  crest  predictions.  Until  this  time,  one 
should  examine  other  alternate  approaches  whicr.  may  cast 
light  on  the  behavior  of  the  response  statistics. 

Although  the  "equivalent  urtimodal  system"  methods  of 
Section  IV  provide  only  an  approx in, at ion  to  the  response 
statistics,  they  appear  to  offer  greater  utility  for  the 
reliability  engineer  than  do  the  more  exact  methods  discussed 
above.  The  concept  of  an  equivalent  >  parameter  makes 
possible  broad  generalizations  concerning  the  response 
statistics.  In  a  qualitative  sense  the  statistics  become 
more  nongaussian  as  v  decreases,  corresponding  to  a  decrease 
in  the  impulsive  rate  or  equivalently  to  a  decrease  in 
system  response  duration  or  an  increase  in  system  bandwidth. 

In  addition  these  methods  yield  a  fair  (order  of  magnitude 
or  better)  quantitative  approximation  for  the  statistics 
on  the  basis  of  very  limited  system  information,  i.e., 

•y/v  and  Thus,  such  order  of  magnitude  predictions  are 
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thought  to  be  invaluable  in  designing  complex  systems  whose 
exact  response  characteristics  cannot  be  knovn  a  priori. 

This  conjectf  t  requires  further  study. 

There  are,  of  course,  some  limitations  to  the  applicability 
of  our  methods  in  the  form  presented  here.  While  the  concept 
of  an  equivalent  system  v  and  the  resulting  qualitative 
statements  are  generally  useful,  for  many  situations  the 
systems  reverberation  time  or  the  bandwidth  defined  by 
Eq .  (41)  are  not  clear.  And,  as  was  mentioned  in  Section  V, 
the  standard  methods  of  measuring  reverberation  time  may 
fail  even  for  a  bimodal  system.  In  this  study  the  rapid 
decay  of  the  impulse  response  of  t  e  bimodal  system  exceeded 
the  writing- speed  limits  of  the  mechanical  chart  recorders. 
Indee.  for  low  Q(Q  <  10)  systems  in  general,  the  utility  of 
a  reverberation  characteristic  is  in  doubt  since  the  decay 
curve  is  not  well-defined.  Rather  than  a  smooth  decay  the 
curve  exhibits  a  series  of  rapid  falls  and  ‘'plateaus"  making 
the  determination  of  a  single  rate  of  decay  rather  arbitrary. 

For  such  systems  some  other  meaningful  measure  of  response 
duration  must  be  developed. 

Like  reverberation  time,  the  bandwidth  concepts  investi¬ 
gated  here  are  limited  in  range  of  applicability.  For 
example,  Eq .  (48),  is  restricted  to  systems  having  highly 
peaked  *-r.-nsfer  characteristics  exhibiting  well-defined, 
re]a»  iy  independent  modes.  For  other  types  of  systems 
the  ir.o-  general  expression  of  Eq.  (40)  would  be  more 
appropriate  since  D  and  T  are  defined  for  all  systems, 
indeed  for  some  systems  an  entirely  different  bandwidth 
definition  might  have  more  utility  especially  if  system 
energies  are  not  readily  available.  In  selecting  a  definition 
of  bandwidth,  however,  one  should  keep  in  mind  that  basically 
it  must  be  a  measure  of  the  inverse  of  character istic  system 
response  time. 

There  also  are  some  difficulties  with  the  particular 
choice  of  ui0,  the  center  frequency  of  a  band  encompassing 
all  significant  modes  of  the  system,  required  to  normalize 
crest  statistics.  For  more  general  systems,  the  choice  of 
such  a  band  of  frequencies  may  not  be  straightforward.  In 
fcuch  casus  perhaps  a  more  satisfactory  choice  for  w0/2v  might 
be  the  zero  -cross ing  rate  for  the  complex  system  excited  by 
Gaussian  white  noise  (e.g,,  Poisson  impulsive  noice  in  the 
limit  v  -*  oo)  .  Thus  the  resonant  frequency  becomes 
(ref.  30)  : 
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with  H(jw)  the  complex  svstem  transfer  characteristic.  Again, 
the  utility  of  definitions  like  Eg.  (50)  must  rely  ultimately 
on  empirical  verification.  It  is  noted  for  the  multimodal 
systems  examined  in  this  studv  that  Eq.  (50)  is  roughly 
equivalent  to  the  "center  frequency"  defined  in  Section  IV. 

Prediction  of  the  appropriate  exact  positive-slope 
zero-crossing  rate,  v+(0),  is  another  problem  that  arises 
when  one  normalizes  the  crest  distribution.  For  overload 
predictions  this  is  unnecessary,  but  for  wearout  prediction 
the  normalized  crest  distribution  is  treated  as  the  "peak" 
distribution.  One  surmises  that  probably  the  damage  over¬ 
estimate  inherent  in  the  use  of  the  normalized  crest  dis¬ 
tribution  (but  the  "peak"  occurrence  frequency  for  the  rises) 
will  persist  even  after  the  extra  factor  wo/2ttv(0)+  is  " 
dropped.  At  any  rate,  more  work  remains  to  be  done  here,  as 
in  the  other  matters  connected  with  the  "equivalent"  unimodal 
model . 


The  entire  development  of  Section  IV  hinqes  on  empirical 
justification.  Thus  while  the  procedures  outlined  for 
selecting  the  "equivalent"  ur.imodal  system  yield  good  agree¬ 
ment  with  multimodal  data,  they  do  not  preclude  the  existence 
of  other,  more  accurate  methods  but  merely  indicate  the 
utility  of  the  concept  of  an  "equivalent"  unimodal  system. 


h-itimittW 


VII 


IMPLICATIONS  FOR  RELIABILITY  PREDICTION 


The  bewildering  diversity  of  failure  modes  encountered 
in  engineering  practice  apparently  can  be  classified  under 
two  headings  (refs.  4,  5):  "overload"  or  "wearout"  failures 

in  terms  of  appropriate  system  response  variables.  Thus, 
in  quite  a  variety  of  cases,  it  is  possible  to  predict  system 
failure  rates  from  mathematical  models  for  the  two  generic 
types  of  failure  mechanisms,  using  knowledge  of  certain 
response  statistics,  through  the  concept  of  a  wear-dependent 
failure  rate. 
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So,  for  overload  failure  occurring  at  response  level 
yc  (A)  ,  which  i:.  general  may  be  a  function  of  total  wear  A 
(at  the  weakest  point),  the  wear-dependent  failure  rate  can 
be  approximated  (refs.  4,  5)  as 

h  { 1 1  A)  =  N+(y  (A)  ;  t)  .  (51) 

o  c 

where  N  ( yc ( A)  ; t)  is  the  average  frequency  of  response  level 
crossings  at  the  weakest  point  in  the  system.  (The  wear- 
dependent  failure  rate,  h(t|A),  can  be  defined  as  the 
probability  density  of  failure  of  thr  system  at  time  t, 
given  perfect  operation  to  date,  and  given  totai  wear  A  at 
said  time.)  Similarly  the  wear-dependent  failure  rate  for 
wearout  failure  can  be  expressed  as  (refs.  4,  5) 

h  ( 1 1  A)  »  v  F  ( A*  -  A)  ,  (52) 

w  1  p 

where  vp  is  the  average  frequency  of  response  peaks,  while 
F(x)  is  the  exceedance  probability  of  partial  wear  (per 
response  "cycle"),  obtained  from  response  "peak"  statistics, 
and  A*  is  the  wearout  level  (usually  unity) .  Generally, 
the  Palmgren-Miner  linear  wear  accumulation  law  is  used, 
and  the  partial  wear,  x,  is  expressed  as  a  power  of  the  peak 
value  S,  i.e., 

x  =  AoSa  (53) 

for  some  A  >0  and  a  >  1,  both  constants  fixed  by  the 
material  (ref.  1).  Alternately,  the  wear  accumulation  procese 
can  be  "linearized"  by  judicious  choice  (ref.  1)  of  expression 
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for  the  analytical  relation  between  "peak"  values  and 
partial  wear. 


To  get  the  overload  and  wearout  failure  rates  proper, 
one  must  average  hQ(t|£)  and  hw(t|A),  using  the  total  wear 
probability  density  for  the  operational  system  in  question. 
Ordinarily,  one  makes  use  of  a  renewal-process-wear  mode.1 
(refs.  32,  33),  which  leads  to  an  asymptotic  normal  wear 
distribution  for  a  system  chosen  at  random,  and  truncated 
normal  wear  (refs.  4,  5)  (truncated  at  £*)  in  an  operational 
system  having  a  definite  wearout  level.  The  failure 
rates  may  be  added  to  obtain  a  conservative  estimate  of  system 
reliability.  Note  that,  all  these  numerical  estimates  using 
wear-dependent  failure  rates — even  apart  from  the  question  of 
finding  more  exact  total  wear  dis tribution--cal 1  for  the 
us  of  a  digital  computer  in  reliability  analysis  (ref.  5). 

To  obviate  the  need  for  extensive  digital  computations, 
one  tries  using  the  "equivalent  single  mode"  approach, 
hoping  to  obtain  "standard"  curves  of  reliability  (wearout 
and  overload  separately,  say)  for  the  impulsive  noise 
environment  in  question  solely  on  the  basis  of  unimodal 
statistics.  In  rough  wearout  reliability  computations, 
of  course,  a  wearout  distribution,  derived  from  the  truncated 
normal  wear  distribution  which  needs  only  knowledge  of  the 
mean  and  mean  square  wear  at  any  time  (in  a  system  chosen 
at  random)  may  be  sufficient.7  But  for  more  refined  estimates, 
without  overescalating  the  computations,  one  would  utilize 
"universal  curves"  or  tables  for  reliability  estimates,  for 
the  "equivalent"  unimodal  system  undergoing  a  statistically 
identical  loading.  Such  data,  presumably,  can  be  readily 
obtained  empirically,  and  thereafter  collected  systematically. 

Unfortunately,  the  situation  is  not  so  cut  and  dried: 
the  anticipated  equivalence  is  not  that  "equivalent"  judging 
by  the  mean  and  mean  square  partial  wear  values  such  as  shown 
in  Table  II  (calculated  for  oy  =  1,  AQ  «=  1,  a  =  4)  .  Here  the 
equivalent  y  and  oj0  values  are  tabulated  for  the  bimodal  and 
multimodal  cases  discussed  earlier  (Figs.  17-9}  .  One  notes 


7It  has  been  found  that  the  averaging  process  introduces  a 
further  conservative  correction  (refs.  4,  5)  in  system  re¬ 
liability  compared  with  using  the  normal  wear  exceedance 
probability  evaluated  at  6*  directly  as  a  wearout  reliability. 
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Table  II  Mean  and  Mean  Square  Partial  Wear 


U> 

o 

<r4> 

8 

<r,  > 

9 

3 

Run  no. 

(if  multi¬ 
modal) 

Y* 

Yd 

2*v+(0) 

0.01 

0.01 

1.61 

28,890 

1 .046xl015 

0.154 

0.367 

119 

0.014 

0.01* 

2.69 

4,215 

8 .696xl010 

0.832 

0.550 

0.05 

0.05 

1.9 

265.3 

2 ,193xl07 

1.13 

0.773 

0.047 

0.075 

1.40 

40.42 

7 

4.283x10 

0.637 

0.689 

118 

0.136 

0.136 

1.05 

58.08 

2 . 939x10"* 

1.16 

0-956 

0.14 

0.17 

1 .46 

19.46 

5 .22xl04 

0.893 

0.909 

117 

0.458 

0.458 

0.57 

12.46 

9.465xl03 

0.807 

1.03 

CD 

00 

1  .00 

1.88 

131.7 

1.414 

2.00 

also  a  fluctuation  in  i  he  normalization  constant  0)^2^  (0) 
for  the  crest  distribution,  indicating  the  aforementioned  em¬ 
pirical  versus  equivalent  unimodal  zero-crossing  rate  discrepancy 
(the  "theoretical"  value  being  .  The  yd  column  indicates 

the  y  yielding  "optimum"  fit,  in  the  experimentalists'  opinion, 
while  y*  is  the  value  obtained  using  the  "bandwidth' procedure 
discussed  in  Section  IV.  The  9  and  &  columns  reveal  that  a 
Weibull  fit  (ref.  5)  was  needed  to  estimate  the  contribution 
of  the  "tail"  (outside  experimental  data  range)  of  the  system 
crest  statistics  to  partial  wear.  Thus,  while  in  the  main  the 
"trends"  seem  correct,  there  are  various  fluctuations  from  a 
"neat"  two-parameter  theory,®  warning  that  the  single  mode 
equivalence  has  limits. 


8  .  + 

Strictly  sj. making,  it  is  not  y  and  o>o,  but  y  and  o>0 /2ttv0  (0) 

that  are  needed  in  wearout  (as  compared  with  overload)  pre¬ 
dictions.  Moreover,  one  needs  an  estimate  for  the  average  fre¬ 
quency  of  peaks  Vp,  for  which  purpose  v+(0)  is  generally  too 
small;  the  average  frequency  of  rises  will  provide  a  conserva¬ 
tive  estimate  for  Vp.  Thus,  the  time  scales  in  any  "standard” 
reliability  curves  must  be  properly  interpreted. 
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To  underscore  the  importance  of  the  "tail",  the  following 
table  gives  the  percent  contribution  of  the  tail  to  the 
multimodal  moments  cited  above;  under  identical  conditions; 


Table  ill  Percent  Contribution  cf  the  Extrapolated 
Region  to  Partial  Wear  Moments 


V* 

Run 

_ l  . 

mean 

mean  square 

0.14 

117 

7.7% 

56.9% 

| 

0.047 

118 

46,8 

9  7.3 

0.01 

119 

99.76 

100.0 

Thus,  the  data  regions  outside  0  <  r,  <  10  increase  in  their 
importance  as  the  response  statistics  become  more  "impulsive" 
in  their  deviation  from  the  Gaussian. 

This  is  further  emphasized  by  comparing  the  results  of  a 
graphi  al  extrapolation  moment  estimation  technique  (utilizing 
a  planimeter)  with  the  Weibull  fit  (corrected  for  deviation 
inside  the  empirical  data  range)  for  the  aforementioned 
cases  (of  Table  II) ;  as  shown  in  Table  IV: 


Table  IV  Wear  Moments:  "Adjusted  Weibull"  versus 
"Graphical  Extrapolation” 


Y* 

O) 

o 

<o4 

> 

<^8> 

2rv+(0) 

AW 

GE 

- - - - 

AW 

1 

GE 

0.01 

1.61 

28, 890 

102  .0 

1 .046xl015 

l.OlxlO5 

0.014 

2.69 

4,  215 

1040 

8.696xl010 

7.7xl06 

0.05 

1.9 

265.3 

224 

2.193xl07 

1.8xl06 

0.047 

1.40 

40.42 

55 . 7 

4.283xl07 

5.29xl04 

0.136 

1.05 

58.08 

49 

2 .939x10  ~ 

4.2xl04 

0.14 

1  .46 

19.46 

36,8 

5.22X104 

2.24xl04 

0.458 

1.00 

12  .42 

12 . 7 

9.465xl03 

7.18xl03 
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The  results  of  the  last  three  tables  speak  for  themselves: 
Clearly,  it  would  be  presumptuous  to  say  too  much  about 
unimodal  equivalence  and  its  validity  for  wearout  (and 
overload)  prediction,  until  the  "tails"  of^the  crest  distri¬ 
butions  have  been  more  thoroughly  studied. 

Various  other  aspects  of  reliability  prediction  might  be 
further  explored  in  this  same  vein,  for  impulsively-driven 
linear  systems,  but  the  above  estimates  of  the  margin  of 
extrapolation  should  drive  home  the  tentativeness  of  these 
considerations:  much  more  work  remains  to  be  done  in  this 

area.  However,  despite  its  limitations,  the  equivalent 
unimodal  approach  seems  to  be  quite  promising  from  an  engineering 
standpoint,  and  until  better  models  for  impulsive  noise  are 
found,  it  can  be  used  for  obtaining  engineering  bounds  on 
system  reliability,  once  a  library  of  unimodal  reliabilities 
has  been  compiled.  (of  course,  some  time  scaling  will 
always  be  necessary,  owing  to  the  ratio  Vp/v+(0)  not  being 
unity,  prior  to  the  multiplication  of  the  wearout  and  overload 
reliabilities  to  yield  system  reliability.)  It  is  hoped 
that  future  investigations  will  further  explore  the  leads 
provided  in  the  present  study,  leading  to  more  complete 
understanding  of  the  perplexing  behavior  of  nongauss ian 
response  processes,  and  their  influence  upon  system  reliability. 


A  two-year  program  for  more  detailed  study  of  unimodal  equiva¬ 
lence  has  been  outlined  by  the  staff  of  the  Propagation  Research 
Laboratory,  Department  of  Electrical  Engineering,  University 
of  Minnesota,  in  its  "Proposal  for  Research  in  System 
Statistical  Response  and  Reliability,"  dated  April  30,  1967, 
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APPENDIX 


CLAMPED  CIRCULAR  PLATE 


There  are  many  possible  systems  which  could  be  used  for 
the  purpose  of  analysis  of  a  multimodal  system  forced  by 
impulse  noise.  The  clamped  circular  plate  was  selected  as 
such  a  system  for  study  for  several  reasons.  FirBt,  it  could 
be  fabricated  with  relative  ease  at  the  machine  shop 
facilities  of  the  Electrical  Engineering  Department, 

University  of  Minnesota.  Second,  the  transfer  characteristics 
of  a  clamped  circular  plate  bear  semblance  to  the  actual 
systems  with  which  this  work  might  be  ultimately  applied,  say, 
fcr  example,  for  fatigue  estimates  of  an  airframe  panel.  Third, 
the  clamped  circular  plate  lends  itself  readily  to  analysis 
(ref.  31) .  This  feature  expedites  the  design  of  the  plate 
and  allows  one  to  place  the  plate  resonant  frequencies  within 
a  workable  range. 

Figure  11  illustrates  the  circular  plate  which  was 
constructed  along  with  its  peripheral  steel  support  ring, 
forcing  assembly,  and  sensing  transducer.  High-leaded  brass 
was  selected  as  the  plate  material  after  an  initial  attempt 
to  fabricate  the  plate  from  aluminum  failed.  The  fundamental 
frequency  was  336  Hz  as  compared  with  the  theoretical  resonance 
of  285  Hz.  The  agreement  between  predicted  resonant  frequency 
and  that  actually  observed  improved  markedly  for  the  higher 
modes . 

The  driving  transducer  was  from  a  Jensen  speaker 
assembly  that  had  been  used  in  previous  unimodal  studies 
(ref.  34)  .  Although  this  forcing  transducer  could  be 
situated  at  four  different  positions,  only  the  one  located 
at  a  distance  of  two  thirds  the  plate  radius  from  the  center 
was  considered  in  the  work  reported  herein.  This  particular 
configuration  seemed  to  excite  the  largest  number  of  modes 
and  therefore  constituted  a  "good"  multimodal  case. 

The  sensing  transducer  is  a  constant-charge  electro¬ 
static  device  which  produces  an  output  voltage  proportional 
to  the  displacement  of  the  plate.  A  Bruel  and  Kjaer  microphone 
amplifier  and  cathode  follower  provides  the  constant-charge 
circuitry  (ref.  35).  The  pickup  probe  (see  Fig.  12),  then, 
can  be  mounted  directly  on  the  end  of  the  B  and  K  cathode 
follower  unit.  A  support  assembly  allows  the  positioning 
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of  the  probe  along  the  center  line  and  rigidly  holds  it 
close  to  the  plate  surface,  say  approximately  0.005  inches 
away.  The  electrical  insulation  isolating  the  support 
assembly  from  the  plate  proper  allows  one  to  determine  via 
a  continuity  check  whether  the  probe  is  positioned 
sufficiently  far  from  the  plate  so  as  to  disallow  physical 
contact  between  the  probe  and  the  plate.  This  electrostatic- 
probe  detection  system  was  found  to  be  more  than  sufficient 
to  detect  accurately  the  forced  vibrations  of  the  plate. 

The  particular  alloy  of  leaded-brass  used  exhibits  very 
little  internal  damping.  This  fact  necessitates  the  applica¬ 
tion  of  damping  tape  to  the  plates  surface  whenever  lower 
modal  quality  factors  are  desired.  The  forcing  transducer 
can  also  provide  some  electromagnetic  damping,  but  this 
effect  is  negligible  when  the  plate  is  forced  at  the  two- 
thirds-  radius  position. 

Table  V  summarizes  the  dimensions  of  the  clamped  circular 
plate  as  well  as  listing  the  material  specifications  for  the 
brass  used.  In  addition,  Table  VI  indicates  the  modes  which 
were  significant  in  the  system  transfer  characteristic  (see 
also  Fig.  6-a-b-c).  The  index  m  designates  the  number  of 
modal  diameters  while  the  number  of  circular  modes  is 
denoted  by  n.  Also  shown  in  Table  VI  are  the  theoretical 
eigenfrequencies  obtained  by  utilizing  the  information  in 
Table  V  . 

Table  V  Physical  Constants  for  Experimental  Plate 


Material:  high- leaded  brass  alloy  #260 

3 

Density:  p  *  8.55  gm/cm 

i 

11  2  1 

Young's  Modulus:  Y  »  11.0  x  10  dynes/cm't  I 

o  i 

I 

Radius:  9.50  inches 

Thickness:  0.092  inches  i 


Table  VI  Resonant  Frequencies  of  Experimental  Plate 


mode  (m,  n) 

0.  1 

■ 

1,  1 

3,  1 

4,  1 

0,  3 

6,  1 

theo.  res  freq.,  Hz 

285 

595 

1420 

1920 

2440 

3210 

actual  res.  freq.,  Hz 

337 

634 

1410 

1890 

2440 

3040 
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Figure  1  Response  of  a  general  time- invariant  linear  system  to  an  impulse 
of  strength  "a"  with  labeling  procedure  for  extrema. 
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Figure  2  Response  crest  distributions  of  a  high  0 
unimodal  system  forced  by  impulsive  noise 
for  several  values  of  the  impulse- 
no  ise-denu i ty  parameter  v. 
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Figure  5  Comparison  of  empirical  response 

crest  distributions  due  to  impulsive 
noise  excitation  of  a  clamped  circular 
plate  under  various  damping  conditions 
with  "generalized-f’  predictions  as 
obtained  by  best  fit. 
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igure  7  comparison  of  empirical  response  crest 
distributions  due  to  impulsive  noise 
forcing  for  the  circular  plate  of 
Fig.  6-b  with  the  approximate  predictions 
based  on  both  system  bandwidth  and 
system  reverberation  time  for  v  ■*  15. 
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Figure  10 


Comparison  of  empirical  response  crest 
dietrj'  *ons  due  to  impulsive  nois< 
forcir,*  tor  a  bimodal  system  with  Qj  =  1C 
Qj  «  20  with  the  approximate  predictions 
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Pigure  11  Illustration  of  the  clamped-circular  plate  construction  indicating 
placesffint  of  the  electromagnetic  driver  and  the  electrostatic 
sensing  probe  mountings. 
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